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Abstract. The authors consider the length, Iff, of the length of the longest 
increasing subsequence of a random permutation of N numbers. The main 
result in this paper is a proof that the distribution function for Iff, suitably 
centered and scaled, converges to the Tracy- Widom distribution [TW1] of the 
largest eigenvalue of a random GUE matrix. The authors also prove con- 
vergence of moments. The proof is based on the steepest decent method for 
Riemann-Hilbert problems, introduced by Dcift and Zhou in 1993 [DZ1] in 
the context of integrable systems. The applicability of the Riemann-Hilbert 
technique depends, in turn, on the determinantal formula of Gessel [Ge] for 
the Poissonization of the distribution function of Iff . 



1. Introduction 

Let Sn be the group of permutations of 1, 2, . . . , N. If it G Sn, we say that 
7r(«i), • • • , 7r(ifc) is an increasing subsequence in it if i\ < %2 < ■ ■ ■ < ik and n(ix) < 
t(*2) < • • • < 7r(ifc). Let In(t) be the length of the longest increasing subsequence. 
For example, if N = 5 and it is the permutation 5 1 3 2 4 (in one-line notation : 
thus tt(1) = 5, 7r(2) = l, ...), then the longest increasing subsequences are 12 4 
and 13 4, and 2jv(tt) = 3. Equip Sn with uniform distribution, 

q n ,N = Prob(l N < n) = -^p, 

where Jn,h = H= (permutations 7r in Sn with In < n). The goal of this paper is to 
determine the asymptotics of q n ,N as N — > oo. This problem was raised by Ulam 
in the early 60's JuJ, and on the basis of Monte Carlo simulations, he conjectured 
that the limit 

(1.1) c= lim -^=E n (In) 

exists. (Here Ejv(-) denotes the expectation value with respect to the distribution 
function q n ,N-) The problem of proving the existence of this limit and the compu- 
tation of c has became known as "Ulam's problem" . An argument of Erdos and 
Szekeres fES[| shows that Ejv(ijv) > — 1, so that if the limit exists, then c > i. 
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Subsequent numerical work by Baer and Brock |BB|] in late 60's suggested that value 
of c is 2. The existence of the limit was rigorously established by Hammersley ]Ha| 
in 1972. In I LS L Logan and Shepp proved that c > 2 and simultaneously Vershik 
and Kerov [VK1| (see also [VK2]) showed that c = 2, thus settling Ulam's prob- 
lem. Alternative proofs of Ulam's problem are due to Aldous and Diaconis [AD|, 
Seppalainen [ |Scl| and Johansson | }Jol| . Over the years, various conjectures have 
been made concerning the variance Uar(Zjv) of In, and Monte Carlo simulations of 
Odlyzko and Rains beginning in 1993, indicated that 



(1.2) 



lim 

N 



Yj^Var{l N ) = c 



for some numerical constant cq ~ 0.819. Also Odlyzko and Rains computed K(In) 
to higher order and found 



(1.3) 



lim 



E(Zjv) - 2VN 



ci 



AT— oo N 1 / 6 

where c\ ~ —1.758. Further historical information on Ulam's problem, together 
with some discussions of the methods used by various authors, can be found in 



[AE| ] and pR | 



Before stating our results, we need to define the Tracy- Widom distribution [ I \\ l| 
(see below). Let u(x) be the solution of the Painleve II (PII) equation, 



OO, 



(1.4) u xx = 2u + xu, and u ~ — Ai(x) 

where Ai is the Airy function. The (global) existence and uniqueness of this solution 
was first established in [HM] : the asymptotics as x — > ±oo are, 

/g-(4/3)s 3 / 2 \ 

u(x) = — Ai(x) + O t-j- I as x — > oo, 

(1.5) 1 



as x 



-oo, 



(see, for example, pMfl , |N[, |DZ2| ). Recall JASJ that Ai(x) 
x — > oo. Define the Tracy- Widom distribution 



2^x 1 / 4 



as 



(1.6) 



F(t) = exp 



(x — t)u 2 (x)dx 



From (jlj) and (JLJ), F'(t) > 0, F(t) 1 as t -> +oo and F(t) -> as t -> -oo, so 
that _F is indeed a distribution function. Our first result concerns the convergence 
of In in distribution after appropriate centering and scaling. 

Theorem 1.1. Let Sn be the group of all permutations of N numbers with uniform 
distribution and let In{^) be the length of the longest increasing subsequence of 
7r 6 Sn- Let x be a random variable whose distribution function is F. Then, as 

N -c oo, 



Xn 



N 



2VN 



X 



i.e. 



lim Prob 



N— >oo 



XN 



'N 



2VN 



ATl/6 



< t 



in distribution, 



F(t) for all t e 
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In order to show that the moments of \n converge to the corresponding moments 
of x we need estimates for the distribution function i*V(i) of xn for large \t\. 
From the large deviation formulas for In (see below), we expect that -Fjv(f) (resp., 
1 — -Fjy(t)) should go to zero rapidly as t — > — oo (resp., t — > +00). In fact, we will 
prove that, for M > sufficiently large, there are positive constants c and C(M) 
such that 



(1.7) 

if -2jV x /3 < i < _ Mj and 
(1.8) 



F N (t) < C(M)e c 



1 - F N (t) < C(M)e 



if M < t < N 5 / e - 27V 1 / 3 . Together with Theorem [□] these estimates yield 
Theorem 1.2. For any m = 1, 2, 3, • • • , we have 

Jim E JV ( X ^)=E( X m ), 

N— +00 

where E(-) denotes expectation with respect to the distribution function F. In par- 
ticular, 

2 



(1.9) 

and 

(1.10) 



lim 

JV— >oo 



Var(Z 



Afl/3 



t 2 dF{t) 



tdF(t) 



lim 



Ejv(Z 



A' / 



2V7V 



ATl/6 



tdF{t). 



If one solves the Painleve II equation (1.4) numerically (see, [ TWl|), and then 
computes the integrals on the RHS of the formulae of ( |1.9| ) and ( |1.10| ), one ob- 
tains the values 0.8132 and -1.7711 which agree with Co and c\ in (1.2) and (1.3) 
respectively, up to two decimal places. 



The distribution function F(t) in Theorems 1.1 and 1.2 first arose in the work 
of Tracy and Widom on the Gaussian Unitary Ensemble (GUE) of random matrix 
theory. In this theory (see, e.g., [Me]), one considers N x N hermitian matrix 
M = (Mij) with probability density 

, N s N 

Z N x e- tr{ - M ^dM = Z N 1 e~ tr ^{^dM il \Y\d{ReM ij )d{ImM ij ), 

where Zn is the normalization constant. In |TWl| , Tracy and Widom showed that 
as the size of the hermitian matrices increases, the distribution of the (properly 
centered and scaled) largest eigenvalue of a random GUE matrix converges pre- 
cisely to F(t) ! In other words, properly centered and scaled, the length of the 
longest increasing subsequence for a permutation ir 6 Sn, behaves statistically for 
large iV like the largest eigenvalue of a random GUE matrix (see the Appendix 
for an intuitive argument). In [TW1], the authors also computed the distribution 
functions of the second, third, • • • largest eigenvalues of such random matrices, and 
the question arises whether such distribution functions describe the statistics of 
quantities identifiable in the random permutation context. 

Recall the Robinson-Schensted correspondence (see, e.g., [ pa| , and also Section 
5.1.4 in |Kn| ) which establishes a bijection it (P(ir),Q(ir)) from Sn to pairs of 
Young tableaux with shape(P(7r)) =shape(Q(7r)). Under this correspondence, the 
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number of boxes in the first row of P(tv) (equivalently Q(k)) is precisely In(tt), 
(see, [65], ]Kn|]). In other words, the results on In can be rephrased as results on 
the statistics of the number of boxes in the first row of Young tableaux. Monte 
Carlo simulations of Odlyzko and Rains |OR[ | indicate that In, the number of boxes 
in the second row of P(tt) (equivalently Q{tt)) behaves statistically for large N, 
like the second largest eigenvalue of a random GUE matrix. More precisely, their 
simulations indicate that 

l N (T N ) - 2y/N 



lira 



and 



lim 



ATl/6 



Var(l N ) 



-3.618, 



= 0.545. 



These values agree, once again, to two decimal places with the mean and variance 
of the suitably centered and scaled second largest eigenvalue of a GUE matrix, as 
computed in [TW1|. Presumably, the number of boxes in the third row of P(ir) 
should behave statistically like the third largest eigenvalue of a GUE matrix as N — > 



oo, etc. In recent work [BDJ|, the authors have shown that this conjecture is indeed 
true for the second row 



Also, beautiful results of Okounkov ]0kJ, using arguments 
from combinatorial topology, have now provided an elegant basis for understanding 
the relationship between the statistics of Young tableaux and the eigenvalues of 
random matrices. Over the last year, many other intriguing results have been 
obtained on a variety of problems arising in mathematics and mathematical physics, 
which are closely related to, or motivated by, the longest increasing subsequence 
problem. We refer the reader to frWj , J|cJ, jjo§ [jo| and pR|. 



As in [ Jo 1 1 , we consider the Poissonization <fi n (X) of q n ,N, 



(1.11) 



0n(A) 



OO 

JV=0 



-Qn.N- 



The function <p n (X) is a distribution function (in n) of a random variable L(X) 
coming from a superadditive process introduced by Hammersley in [Ha|, and used 
by him to show that the limit (1.1) exists. The random variable L(X) is defined 
as follows. Consider a homogeneous rate one Poisson process in the plane and 
let L(X) denote the maximum number of points in an up-right (increasing) path 
through the points starting at (0, 0) and ending at (vA, V~X). For more details see 



1.1 



|ADj] and [3c2|, and for a generalization to the non-homogeneous case see [DeZel 
Theorem 



and 1.2 hold for the random variable L(X) as A — > oo. Referring to the 
"de-Poissonization" Lemmas 3.2 and 8.3 below, we see that it is easy to recover the 
asymptotics of q n ,N as N — > oo from the knowledge of <f> n {X) for A ~ N. In other 
words, in order to compute the asymptotics of In, we must investigate the double 
scaling limit of 4> n {X) when A — * oo and 1 < n < N ~ A, and this is the technical 
thrust of the paper. 

To this end we use the following representation for <f> n {X), 

(1.12) (f> n (X) = e- A J D„-i(exp(2VAcos0)), 

where £>„_i denotes the n x n Toeplitz determinant with weight function f(e t6 ) = 
exp(2\/A cos 8) on the unit circle, (see, e.g.,[Szl|). The above formula follows from 
work of Gessel in [Ge| using well known results about Toeplitz determinants. As 
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noted in [ Jolj, the formula can also be proved using the following representation 
for q niN , l<n<N, discovered by pPWW| , 



(1.13) 



Qn,N 



2 2N N\ 
(27V)! 



(5^COS0j 



,2N 



n \^ 

l<j<k<n 



(2n) n n\ 



In addition, an earlier result of Diaconis and Shahshahani ([DS]) shows that the 
above formula ( 1.13| ) is true also in the case n > N when q n ^ = 1. Inserting (1.13) 
into (1.11), we obtain 



(1.14) 4> n {\) =e 



— A 



1 r " 

^/ e xp( 2VAEcos^) n 



^ififc 1 2 m, 



d n 9, 



which is p recisely ( |l,12|) by standard m etho ds in the theory of Tocplitz determinants 
(see, pzjj ). An additional proof of (U.12 ) can be found in [ G WW | , and also an 
alternative derivation of formula ( 1.13| ) is given in |Ra]. For the convenience of the 
reader we provide (yet another) proof of (1.12) in the Appendix to this paper. 

Using the integral representation ( 1.12j ), Johansson ( [[lolf ) proved the following 
bound for (f>(\) : for any given e > 0, there exist C and 8 > such that 



(1.15) 



< <j> n {\) < Ce~ sx if (1 + e)n < 2V% 
C 

O<1-0„(A)<— if (l-e)n>2V\. 



This information and the de-Poissonization Lemma S.2 are enough to give a new 
proof (poll) that 



(1.16) 



lim L N /2VN = 1. 

AT— *oc 



The first estimate in (1.15) is a consequence of the following lower tail large devia- 
tion formula for </>„(A), 



(1.17) 



1 



lim — d>, /ti (A) = — 1 + 2x - 



"—X 2 — ^— log — 



-U(x), 



if x < 2. For the upper tail Seppalalnen i n ]Se2j used the interacting particle system 
implicitly introduced by Hammersley in fHafl to show that 



(1.18) 



lim log(l - (A)) = -2xcosh" 1 (j;/2) + 2\/x 2 - 4 = -I{x) 



if x > 2. We note that Hammer sley 's interacting particle system was also used ear- 
lier by Aldous and Diaconis in |AD|] . The sup er-additivity of the process described 
above implies that we actually have, see Se2] and also []kJ, 



(1.19) 



l-4> [xM] (M 2 )<e- MI ^\ 



if M is a positive integer and x > 2. This estimate can be used to show ( 
but in this paper we will give an independent proof of ( fLS| ). The large deviation 
formula (1.18) implies, via a de-Poissonization argument, that for x > 2, 

(1.20) 



lim — — \ogProb(l N > xVN) = —I(x) 
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For the lower tail the large deviation formula for In is not the same as for L(X), 
the Poissonized case. Deuschel and Zeitouni in [DeZeS] use combinatorial and 
variational ideas from Logan and Shepp JLS| to prove that 



(1.21) 
if < x < 2, where 

(1.22) H(x) 



N 



lim — log Prob{l N < xVn) 



-H(x), 



2x 2 



1 X X X 

'2 + T + log 2-( 1 + T) log ( 4 + ^ 



For the lower tail we have no analogue of fll.lSp . The rate functions U and H are 
related via a Legendre transform, see |Se2]. The above results show clearly that 
the distribution function for In is sharply concentrated in the region {(2 — e)^fN < 
In < (2 + e)VN} for any e > 0, and they can be used to see heuristically that the 



variance for In should be of order iV 1 / 3 , see [Ki 



As is well known (see, Szl]) the Toeplitz determinant D n -i in ( 1.12| ) is inti- 
mately connected with the polynomials p n (z;X) — n n (X)z n + which are or- 
thonormal with respect to the weight f(e l6 )^ = exp(-\/A(z + z -1 )) ^j- on the unit 
circle, 



(1.23) 



(W 



Pn{e l »)p m {e ie )f{e w )— = S n , 

Z7T 



for n, m > 0. 



The leading coefficient fc„(A) can be expressed in terms of Toeplitz determinants, 



(1.24) 



<(A) 



D n (X) 



where D n (\) = D n (exp(2^/\cos9)). But by Szego's strong limit theorem ([ 3z2|) 
for Toeplitz determinants, lim n ^oo D n (X) — e A , and hence 



(1.25) 



log^ n (A) = ^log4(A). 



Therefore, if one can control the large k, A behavior of k|(A) for all k > n, one will 
control the large n, X behavior of <j> n (X). 

The key point in our analysis is that k\(X) can be expressed in terms of the 
following Riemann-Hilbert Problem (RHP) : Let E be the unit circle oriented coun- 
terclockwise. Let Y(z; k + 1, A) be the 2x2 matrix- valued function satisfying 



(1.26) 



'Y(z; k + 1, A) is analytic in 
Y+(z;k + l,X)=Y-(z;k+l,X) 



_i_„VA(z+z-y 

7?+T e 



- s 

v° 1 

Y(z;k + l,A)z-( fc+1 ) <T3 = / + 0( i ) as z -> oo, 



on S, 



where Y + and F_ denote the limit from inside and outside of the circle respectively, 



and (73 = ( o _i ) , so that z 



-(k+l)a 3 



Here / is the 2x2 identity 



matrix. This RHP has a unique solution (see (4T) below), and the fact of the 
matter is that 



(1.27) 



4(X) =-F 2 i(0;fc + l,A) 
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where >2i(0;fc, A) is the (21)-entry of the solution Y at z = 0. In [DZ1| and 
|DZ2 ], Deift and Zhou introduced a steepest descent type method to compute the 
asymptotic behavior of RHP's containing large oscillatory and/or exponentially 
growing/decaying factors as in (1.26). This method was further extended in [ DVZl| 
and eventually placed in a very general form by Deift, Zhou and Venakides in 
DVZ2fl , making possible the analysis of the limiting behavior of a large variety of 



asymptotic problems in pure and applied mathematics (see, e.g., |DIZ|). As we will 
see, the application of this method to (1.26) makes it possible to control the large 



k, A behavior of k£(A). The calculation in this paper have many similarities to the 
computations in [DKMVZ], where the authors use the steepest descent method to 
obtain Plancherel-Rotach type asymptotics for polynomials orthogonal with respect 
of varying weights, e~ NV ^dx on the real line, and hence to prove universality for 
a class of random matrix models. The Riemann-Hilbert formulation o f the theory 
of ort hogon al polynomials on the line is due to Fok as, Its and Kitaev ([ FIK |) : the 
RHP ( 1.26| ) is an adaptation of the construction in [ FIK to the case of orthogonal 
polynomial with respect to a weight on the unit circle. 

This paper is arranged as follows. In Section 0, we discuss of some of the basic 
theory of RHP's and also provide some information on the RHP associated with 
the PII equation. This information will be used in t he construction of an approx- 
imate solution, i.e. a parametrix, for the RHP ( 1.26 ) in subsequent sections. The 
appearance of the PII equation in the limiting distribution F(t) for xn originates 
in this construction of the parametrix. A connection of </> n (A) to Toda lattice and 
the Painleve III equation is p resen ted in Section [| S ectio n || is the starting point 
for the analysis of the RHP ( 1.26 ). In this section, ( |l.26| ) is transformed into an 
equivalent RHP via a so-called (/-function. The role of ^-function, first introduced 
in |pZ2)| , and then analyzed in full generality in [ DVZ2 1 , is to replace exponentially 
growing terms in a RHP by oscillatory or exponentially decreasing terms. It turns 
out that in the case of ( 1.26 ), as in DKMVZ 1 , the g- function can be constructed 
in terms of an associated equilibrium measure dfi(s) as follows, 



(1.28) 



9(z) 



log(z — s)dfx(s) 



The measure d/j, is the unique minimizer of the following variational problem 



(1.29) 
where 
(1.30) 



E v = inf{/ v (fj,) : /i is a probability measure on the unit circle E} 



log|s — w\ 1 dfi(s)dfi(w) 



ExS 



V(s)djl(s) 



The variational problem (1.29) describes the equilibrium 



and V(s) = — \/A(s+s~ 
configuration of electrons, say, confined to the unit circle with Coulomb interactions, 
and acted on by an external field V. It turns out that the support of the equilibrium 
measure depends critically on the quantity 

We need to distinguish these two cases, 7 < 1 and 7 > 1. As noted by Gross and 
Witten (| GW |), and also by Johansson ([ [Jol| ), the point j^j = 1 corresponds to a 
(third order) phase transition for a statistical system with partition function (1.14). 
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The first case, when 7 = 



< 1, is discussed in Section [| and the second case, 

when 7 = > 1, is discussed in Secti on ^ . The principal results of the above 
two sections are summarized in Lemmas t5.ll and |6.3|. We obtain full asymptotics 



of k„(A) for n, A > when n, A — > 00. In Section 
k > n, we obtain the asymptotics of </> n (A) in Lemma 7.1 
4> n {N) and g„jy (de-Poissonization) is discussed in Section ||. Finally, the proofs 



by summing up k|(A) for all 
The relation between 



of Theorem 1.1 and 1.2 are given in Section 



Notational remarks : The primary variables in this paper are n, N, and A. 
The letters C, c denote general positive constants. Rather than introducing many 
such constants Ci, C2, • • • , c\, C2, • • • , we always interpret C, c in a general way. For 
example, we write \f(x)\ < 2C|.g(a;)|+e c |/i(a;)| < C(\g(x)\ + \h(x)\), etc. We will also 
use certain auxiliary positive parameters M, Mi, M2, • • • , M7. If a constant depends 
on some of these parameters, we indicate this explicitly, for example, C(M2,M&). 
In addition to the standard big O notation, we also use a notation Om- Thus 
/ = 0(^x73) means |/| < t^t/J: where C is independent of M, Mi, ■ ■ ■ . On the other 



hand, / = M {-z^m) means |/| < 



C(M,Mi,--- ) 



depends on at 



7T73 '-, where C{M,M X , 

least one of the parameters M, Mi , ■ ■ ■ . 

In the estimates that follow we will often claim that an inequality is true 
n — > 00". For example, in (7.3) below, we say that 

2VA 



as 



I log K (A) I < Cexp^-c(n + l)(l- 



,3/2 



V 



as n — > 00. This mean that there exists a number no, say, which may depend on 
all the other relevant constants in the problem, such that the inequality is true for 
fi > no, etc. (For this particular inequality the only other parameter is M 5 , but it 
turns out that the constants C, c can be chosen independent of M5 (see below).) 

2. Riemann-Hilbert Theory 

In this section, we first summarize some basic facts about RHP's in general, and 
then discus s th e RHP for the PII equation in some detail. Basic references for 
RHP's are |x|, |GK| , and the material on PII is taken from |DZ2 ] . 

Let E be an oriented curve in the plane (see, for example, Figure By con- 




Figure 1. 



vention, the (+)-side (resp., (— )-side) of an arc in E lies to the left (resp., right) 
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as one traverses the arc in the direction of the orientation. Thus, corresponding to 
Figure [l], we have Figure ||. Let Eq = E — {points of self-intersection} and v be a 




Figure 2. 



smooth map from So — > Gl(n, C), for some n. If E is unbounded, we require that 
v(z) — ► I as z — > oo along E. The RHP (£, v) consists of the following (see, e.g., 
CG j) : establish the existence and uniqueness of an n x n matrix valued function 
Y(z) (the solution of the RHP (£, v)) such that 



(2.1) 



Y(z) is analytic C — E, 
Y+(z) =Y-(z)v(z) zeE , 
Y(z) —¥ I as z — > oo. 



Here i±(z) = liru Z '-»2 where z' G (±) — side of E. The precise sense in which 

these boundary values are attained, and also the precise sense in which Y(z) — > / 
as z — » oo, are technical matters that should be specified for any given RHP (E, v). 
In this paper, by a solution Y of a RHP (E, v), we always mean that 

Y(z) is analytic in C — E and continuous up to the boundary 
(including the points in E — Eo) in each component. 
The jump relation Y + (z) = Y-(z)v(z) is taken in the sense of 
continuous boundary values, and Y(z) — * I as z — > oo means 

Y(z) = I + 0(-j — r) uniformly as z — > oo in C — E. 



(2.2) 



Given (E, v), the existence of Y under appropriate technical assumptions on E and 
v, is in general a subtle and difficult questio n. Ho wever, for the RHP ( 1.26|) , and 
hence for all RHP's obtained by deforming ( 1.26| ) (see, e.g., ([Q)|)), we will prove 
the existence of Y directly by construction (see, Lemma 4.1) : uniqueness, as we 
will see, is a simple matter. 

The solution of a RHP (S, v) can be expressed in terms of the solution of an 
associated singular integral equation on E (see, (|]t]), ( |2.8| ) below) as follows. Let 
C± be the Cauchy operators 



(2.3) 



(C±f)(z)= lim 



f(s) ds 
s — z' 2iti ' 



z G E, 
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where z' — > z± denotes the non-tangential limit from the ± side of S respectively. 
A useful reference for Cauchy operators on curves which may have points of self- 
intersection is [GK|. Under mild assumptions on E, which will always be satisfied 
for the curves that arise in this paper, the non-tangential limits in ( |2.3| ) will exist 
pointwise a.e. on £. Furthermore, if / G L P (T,, \dz\), 1 < p < oo, then the bound- 
ary values (appropriately interpreted at the points £ — £o of self-intersection) of 

Je T=l4k are als0 taken in the SenSe ° f LP aild \\ C ±f\\Lv(Y,,\dz\) < C p \\f\\ L P(j:,\dz\)- 

A simple calculation shows that 

(2.4) C+-C- = 1. 
Let 

(2.5) v = bZ 1 b + = (I -w-)~ x (l + w+) 

be any factorization of v. We assume 6±, and hence w±, are smooth on So, and if 
E is unbounded, we assume b±{z) — > I as z — > oo along £. Define the operator 

(2.6) C w (f) = C+(fw-) + C-ifw+). 

By the above discussion, if w± € L°°(E, \dz\), then C w is bounded from L 2 (E, \dz\) — ► 
L 2 (E, Suppose that the equation 

(2.7) (1 - C w )fi = 1 on E 

has a solution fi 6 7 + £ 2 (E), Or more precisely, suppose ^ — I € ^ 2 (^) solves 

(2.8) (1 - C w )(fx -I) = C W I = C+{w^) + C-(w+), 

which is a well-defined equatio n in £ 2 (E) provided that w± £ L°° n L 2 (T,, \dz\). 
Then the solution of the RHP (fy]) is given by (see, |CG|,|bc|) 

(2.9) Y{z) =I+[ MM') z i s . 



Indeed for a.e. z € E, from (2/7) and ( |2.4| ), 

F+(z) = J + C + ( M (s)(^+(s) +te-(a))) 
= I + C w (p) + (C+ -C-)(jiw+) 
= fi + fiw+ 
= v{z)b+(z), 

and similarly Y_(z) = /i(z)6_(z), so that l+(z) = y_ (z)6l 1 (z)6 + (z) = F_(z)u(z) 
a.e. on E. Under the appropriate regularity assumpti ons on E and v, one then 
shows that Y(z) solves the RHP (E, v) in the sense of (p"^). 

As indicated, the above approach to the RHP goes through for any factorization 
v = (I — w_) _1 (J + w+). Different factorization may be used at different points in 
the analysis of any given problem (see e.g. [DZ1]). However, in this paper we will 
always take W- — 0, so that v = (I + w + ). Thus C w always denotes the operator 
C-(-(v-I)). 

In this paper we will not develop the general theory for the solution of RHP's, 
giving conditions under which (2.7) has a (unique) solution, etc. Rather, for the 
convenience of the reader who may not be familiar with Riemann-Hilbert theory, 
we will use the above calculations and computations as a guide, and verify all the 
steps directly as they arise. 
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We now consider the RHP for the PII equation (|FN|], |JMU[ : see also |n|, 



FZ], |DZ2|). We will consider two equivalent versions of the RHP for PII. These 



two RHP's will be used in the later sections for the construction of parametrices 



for the solution of ( L26| ). 

Let £ P7/ denote the oriented contour consisting of 6 rays in Figure |. Thus 
£ p// = U^ =1 {S P// = e^-^^R+j, with associated jump matrix v PH : £ p// -> 
A7 2 (C), where the monodromy data p, q and r are complex numbers satisfying the 
relation 



(2.10) 



p + q + r + pqr = 0. 



For x € R and z € S p// - {0}, set 




Figure 3. v pn and T, pu 



(2.11) 

where 
(2.12) 



,pii. 



i ■ . (z) = e ™ r "-'v' "e" 



4z 6 



ypn 



+ xz. 



0. 



The contour S P7/ consists precisely of the set i?e(i4z 3 /3) 
particular, that v PII (z) — I ^ L 2 (E P// ). For example, as z 
axis, w P// (z) — / is oscillatory (on the other rays, S P// , k = 
could grow), and so we cannot expect that the RHP (E p// , w P// ) has a solution 
in the sense of (2.2). However, if we rotate E P7/ in the clockwise direction by any 
angle 6 , < 6 < tt/3, S p// S p// = e- i9 °T, PI1 ', then it is easy to see that 



This implies, in 
-> +oo along the real 
1,2,4,5, v PII (z)-I 



,piii 



(z) - 7 e L 2 n L°°(S P77 ), and we may expect that the RHP (S PJ7 , v pu ) has 



a solution in the sense of ( |2.2| ). Moreover, as v p (z) is analytic, it is clear that if 
one can solve (T,^ 1 ,v PI1 ) for some < 8 < tt/3, then one can solve (E PI1 ,v Pn ) 

for any other < 6*o < 7r/3, and the solution of the 6*o-problem can be obtained 
from the #o-P r oblem by an analytic continuation, and vice versa. So suppose that 



12 



JINHO BAIK, PERCY DEIFT, AND KURT JOHANSSON 



for some fixed < 6q < tt/3, and for mg f ^ I (z; x) is a (2 x 2 matrix) solution 

of the RHP (Sf o 7/ , // ), 



- £ 



(2.13) 



m,g (z) analytic in 

+ = W)J*)«f"(*). ^ z G 
nig 11 (z) —* I as z — > oo, 



in the sense of (2.2). Let mf II (x) denote the residue at oo of m,g (z), given by 



<-(z;,)=/ + <^M + 0(l) 

as z — > oo. Then 

(2.14) u(x) = 2im^[ I 2 (x) = -2im" 2 [(x) 

solves PII (see |]|], pMU| ), 



u xx = 2u + xu, x £ 



where mf [ 2 (x) (resp., mf^Kx)) denotes the (12)-entry (resp, (21)-entry) ofmf 



phi 



It is easy to see that mf 7/ (a;), and hence u{x) in ( [2.14 ), is independent of the choice 
of 6»o € (0, tt/3) 



A solution of the RHP (E 



pii ..pii 

So ' % 



) for some 0q, hence for all 9q € (0, 7r/3), may 



not exist for all p, q, r satisfying ( 2.1C ) and ieK. A sufficient condition (see |FZ|) 
for the RHP to have a unique solution (in the sense of (^^)) for all x E R, is that 

|<?-p| < 2 and rel. 

In this paper, we need the singular case 

(2.15) p=-q = l and r = 0. 

The latter condition r = implies that there is no jump across the rays ± e l ( 27r / 3 ~ o) j 




FIGURE 4. p = -g = l,r = case ; E^ 77 ' 1 and v PI1 ' 1 



and we may replace E^, 77 by E 7 ^ 77 ' 1 as in Figure || (note that the orientations across 
the rays e~ l6 ° , e l ( 2,r / 3 - 8 o) have been reversed). As noted in DZ2], a unique solution 
in the sense of (2.2) still exists in this singular case for all x € R : a proof of this 
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fact is not given in [DZ2], but can be found in [DKMVZ3 ; nonregular case, Case 
II] . In addition, the solution has the property that 



m^ 77,1 (z; x) and its inverse are uniformly bounded 

for (z, x) G (C • 
pii,i, 



(2.16) 



E™)x[-M,M], 



for any fixed M > 0. As m^ ) /7 ' 1 (z; x) solves ( [2.13| ) in the sense of ( |2.2| ), we see in 
particular that ( 2.16 ) holds up to the boundary in each sector. 



The asymptotics of u(x) = 2im P { 7 ! ' 1 {x) given in ( |l.5| ), is computed in [ DZ2 | via 
the above RHP and from the proof in |DZ2| , one learns that 



pii,i 



(2.17) 



where 



(x)=0 



-(4/3):r 3/2 



pii.i, \ 

m 1.22 W 



- x as x 



as x — ► oo, 



-oo, 



pii.i 



denotes the (22)-entry of mf . Also, using the methods in |DZ2|, 



for example, one obtains the relation 



(2.18) 



2im 1 j' I 2 I 2 ' 1 (x) = u 2 (x). 



dx 



and verifies directly that 2imf 22 ,1 ( a; ) is real- valued 




Figure 5. £ P7p2 and v PII > 2 

For the first of the two equivalent RHP's advertised above, we consider Figure |[ 
which consists of the real axis (the dotted line), S^ 77 ' 1 for some fixed, small 6q > 
(the dashed lines), and a contour E P77,2 consisting of a pair of curved solid 
lines. The contour E P77,2 is of the general shape indicated in the Figure, with 
one component in C+ and one component in C_, and we require that E P77,2 is 
asymptotic to straight lines lying strictly within the region {|ar<7z| < 7r/3}U{27r/3 < 
argz < An/ '3}. Together with the line {xe~ l6 ° : x € R}, these contours divide the 
complex plane into 4 open regions, fl k ' 2 , k = 1, 2, 3, 4, as shown in Figure 0. Let 
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PII.2 



be the jump matrix on £ p/p2 which is given by ( \ \ ) in C+ and by ( J j 1 ) 



v 

in C_ . We define 



m 



m 



PII,2 



,PII.2 



,P//,2 



.P//.2 



m PII,l e . l{ g Pll)ada3 ^oyl . 
m W/,l e -,( fipn )„ dj3 ji„j 
PII,l p -i{e PU )ada 3 ^l -I) -1 



m 



m 1 2 X PI sZ 2 j 

r,PII.2 „ oP^.l 
m S2 2 (I i 2 X , 

m S 2 3 n s 2 4 , 
r,P//.2 _ P^,1 



TO fl otherwise, 



where the regions ri P7pl ,/c — 1,2,3,4 are defined in Figure |4|. A straightforward 
calculation with the jump relations for mf 11 ' 1 , shows that m PI1,2 solves the new 
RHP 



(2.19) 



to p// > 2 is analytic in C - £ p// > 2 , 

PiT.2 P77,2 pa 2 ^PII 2 

to , = to_ ' v~ ' on L ' , 



,P/P2 



/ as z 



where f P77 ' 2 = e l { e Pu) ad(7 3 V Pii,2 anc j v pii,2 j s gj ven m Figured. This deformed 

P II 1 

RHP is clearly equivalent to the original RHP for m 6o ' in the sense that a solution 
of the one RHP yields a solution of the other RHP, and vice versa. Also we have 



(2.20) 



TO-, 



PII.2 



for the residues of nig 11 ' 1 (resp, m PIP2 ) at oo. From (2.16), we see that for any 
fixed M > 0, 



(2.21) 



TO 



P//.2 



(z; n) and its inverse are uniformly bounded 



for (z, x)e(C- £ P7p2 ) x [-M, M]. 
A particular choice of contour £ p// will be made in Section || (see below) . 




Figure 6. S p/p3 

The second of the equivalent RHP's is restricted to the case x < 0, and we 
consider Figure |[ which consists of the real axis (the dotted line), S^ 77 ' 1 for some 
fixed small 9 > (the dashed lines) and a contour £ P77 ' 3 = uf. =1 E P77 ' 3 consisting 
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PI I 3 

of 5 straight lines, one finite and four infinite. The regions fl k ' , 1 < A; < 4, are 
the components of C — E p// ' 3 . 

The infinite lines make an angle strictly between and 7r/3 with the real axis. 
Set 

3 2 



(2-22) ff PW W = i(* a + ^r 



Y, J -^], and behaves like |z 3 + xz 



which is defined to be analytic in C — [ 

|1 + O(js) = 0pii(z) + 0{\) asz^oo. Therefore for any M > 



(2.23) 

e i(g pn {z)-0 P „(z)) 

and 

(2.24) e" 
We define m P7p3 by 



is bounded for (z,x) € (C — [— 



2 ' 



])x[-M,0] 



1 as z — > oo uniformly for — M < a; < 0. 



m P//,1 [e~ i(ep " )a<icr3 ( } 5)] e i (f P "( 2 )- e P"(z))o'3 

P",lf p -i(epji)o<fcr3(l -1 y 1 ] p i(g PII (z)~9 PII (z))a 3 



r,PII,3 r-.PII,l 

m ii 3 n sZj , 



'On 



m 



■' ! r _-i(e P „)adcr3 ( 1 -1 )] e i(s P "(z)-ep// (z))o- 3 



'00 

m e ' otherwise 



in il 
in il 



F "' 3 nnf"-' 
n(n£"-' u<>' : 



2 4 

P7/,3 „ mPIl.l , , ,-,P77,^ 



Then from the jump relations for wig* ' 1 , we see that m PI1,3 solves the new RHP 



(E P7p3 , W P77 - 3 ) in the sense of (gj), 



,P/i,3 



is analytic in C — E 



PII,3 



(2.25) 



P77,3 P77,3 P77 3 vP/J 3 

I , = 771 _ Vt. '° on E ' , 



,P77,3 



/ as z — > oo, 



where v PI1 ' 3 is given by 



(2.26) 

Also we have 
(2.27) 



1 



1 -e~ 2 ^ 
,0 1 

1 

PII 
m n = m-. 



v.P/7,3 ^P77,3 

on ,E 2 



^P/7.3 ^P77,3 

on E 3 ,E 4 



on E 



PIP3 



P77,3 



0"3, 



for the respective residues of m e ' and m PIP3 at oo. Finally, from ( 2.16] ) and ( 2.23| ), 
we see that, for any fixed Mel 

m PII ' 3 (z; x) and its inverse are uniformly bounded 

^ 2 ' 28 ^ for (z, x) 6 (C - E p/P3 ) x [-M, 0] . 
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3. Connection to the Toda Lattice and the Painleve III Equation 



In this Section, we discuss the connection of the RHP ( 1.26[ ) for k 2 and the 
RHP for the Toda lattice and the Painleve III equation. In the RH context, the 
connection results from the specific form of the weight, e ^( z + z ) . Connections can 
also be seen from the Toeplitz determinant /orthogonal polynomial point of view as 
in |P|], |H and fw|. The purpose of this short Section is purely to establish the 
various connection, but we do not use the results in the sequel. 



Write q = k + 1 in ( |1.26[ ). We define 

)Y(z;q)\ 



(3.1) m TL (z;q) 








o 











\z\ < 1. 



1 
-1 



1*1 > 1, 



)Y(z-q) 

A simple calculation shows that m TL solves the following RHP. 
m TL (z) is analytic in C — E, 

ml L {z)=m T _ L {z) \^ zqeV i {z _ z -i } 

m TL (z) —> I as z 



(3.2) 



z -q e -V\(z- 
1 



Once again, the RHP for Y is equivalent to the RHP for m TL in the sense that a 
solution of one problem yields a solution of the other problem. 

Recall that the RHP related to the Toda Lattice problem, for — oo < m < oo, 

^=2(6^-6^) 

(3.3) * 

" ,u m i t \ 

— 7— — «m 8ra+l — flmji 

at 



under initial data decaying at infinity is the following (see,e.g.,[Ka]). Suppose that 
there are no solitons and denote the reflection coefficient by r(z), z £ E. Then we 
find Q(z) such that 



(3.4) 




E 

l-| r (z)|2 _ f ^ z 2m e -t( Z -z-^ 

r(z)z- 2m e t( - z - z ~^ 1 
: — > oo. 



When q is even, if we set y/\ = t and q = —2m in (3.2), then the RHP is identical 
with the above RHP with r(z) = 1. 

For the connection to the Painleve III equation, define 



(3.5) 



phi 



'(-l)% Ti (z) 
m TL (z) 



\z\<l, 
\z\ > 1. 



Note in (gj), 
(3.6) 



(-1)* 







_ 7 -9„-v / A(z-z 1 ) 



z H e y 



1 IN 

-1 0; 
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where z' is analytic in C — (—00, 0] and real- valued for real z. If we set \[\ - 



-ix, 



then this is the same RHP for the particular Painleve III equation (see [FMZ] for 
results and notations) 



U 



-U X + -{-Aqu 2 + 4(1 - qj) 

X X 



* 3 - 4 
Au -\ 

u 



with monodromy data 



Voo — —Uo — 9) 

»0 = b() = doc = 

E 



0. 



-Vo)- 



In the RHP (1.26), we are interested directly in the quantity — l2i(0; k + 1, A), 
or by (3.1), mf^O; q). On the other hand, for the Toda lattice and the PHI 
equation, one is interested in quantities other than m|] L (0;g) which are related to 

the respective RHP's. For example, the solution u(x) of PHI equation is given by 

pi 1 1 

u(x) = — ixirrii III )i2 where m PIH = I + h 0{\), which is clearly different 

from (— l) q m^ n (0; q). However, the importance of the connection of ( 1.26| ) to the 
RHP's for Toda lattice and the PHI equation lies precisely in the fact that (a m , b m ) 
(resp., u(x)) solve differential-difference (resp., differential) equations which in turn 
imply that the coefficients of the generating function </>„(A), 2a/A = — ix, must 
satisfy a certain class of identities. We plan to investigate these relations in a later 
publication. 

Finally, note that for PHI, the interesting asymptotic question is to evaluate the 
limit x = — > 00, with q fixed. In this paper, as in the Toda lattice, we are 
interested in the double limit when A — > 00 and q is allowed to vary (note that in 
Ka], the singular case r(z) = 1 is not considered). When A — * 00, ~ 1, we 



are in a region where the solution of PHI equation degenerates to a solution of PII 
equation, and this explains the appearance of PII in the parametrix for the solution 



of Y of the RHP (1.26) 



4. Equilibrium Measure and g-function 

In thi s Se ction, the equilibrium measure is explicitly calculat ed f or each 7 > 
(Lemma O) and, using this e quilib rium measure, the (/-function (|4.8| ) is introduced 
in order to convert the RHP (|L26b into a RHP which is normalized to be / at 00. 



Let E denote the unit circle oriented counterclockwise and f(e l °) = f(z) be a 
non-negative, periodic, smooth function on S. Let p q {z) = n q z q + • • • be the q- 
th normalized orthogonal polynomial with respect to the weight f(e te )^ on the 
unit circle. Define the polynomial p* q {z) = z q p q (l/z) = z q p q (\/z) (see [ |5zl| ). We 
consider the following RHP : Let Y(z) be the 2x2 matrix- valued function satisfying 



(4.1) 



is analytic in C — E, 

Y(z)z-'- k+1 ^ 3 =I + 0{\) as z^oo. 
The following Lemma is the starting point of our calculations. 
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Lemma 4.1. (cf. jFIKj , |DKMVZl| ]j The RHP (Q has a unique solution 



Proof. Existence : Using the property of Cauchy operator C+ — C_ = /, where 
Ch{z) = J s ^jt~s, it is a straightforward calculation to show that the above ex- 
pression for Y satisfies the jump condition. The asymptotics at oo codes in precisely 
the fact that the p' k s are the normalized orthogonal polynomials for the weight 
f(e l9 )d8. 

Uniqueness : Suppose that there is another solution Y of RHP. Noting 

det i ^/« = i 

, we have that det Y is entire, and — > 1 as z — > oo. Therefore 

v° 1 / 

by Liouville's theorem, detF =1. In particular, Y is invertible. Now set Z = 
YY^ 1 . Then it has no jump on S hence is entire. Also, Z — > I as z — > oo, and 
therefore Z = I. □ 

From this Lemma, we have 

(4.2) K 2 = _ y2i(Q) 



Therefore the RHP ( 1.26 ) has a unique solution and ( 1.27 ) is verified 



Again set q = k + 1 and 
(4.3) T -2£ 

We are interested in the case when q and 2-\/A are of the same order , or more 
precis ely, 7 — * 1. In this Section, and also in Sections || and S, we consider the 



RHP ( 1.26 ) with parameter 7 and q, 

{Y(z; q) analytic in C — S, 
Y+(z;q) 
Y(z;q) = (I + 0(i))z« ff3 as z -> 00, 

rather than A and g. With 7 fixed, the RHP (|P| ) is of the Plancherel-Rotach 
type with varying exponential weight e 3 r^ z+z ) on the unit circle (see pzl| , 
DKMVZ~L| ) . Similar problem on the real line is analyzed in [DKMVZ1] with- 




out double scaling limit (7 is kept fixed). Our goal is to find the large q behavior 
of y 2 i(0;g) for all 7 > 0. 




Figure 7. branch cut of log (2 — e ie ) 
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Let dfx(s) be a probability measure on the unit circle. Define 
(4.5) g(z)= / log(> - s)d(j.(s) 



where for each 9, the branch is chosen such that log(z — e l6 ) is analytic in C — 
(— oo, — 1] U {e lt : — tt < t < 9} (see Figure [?]) and log(z — e 10 ) ~ logz for real 
z — > oo . The following Lemma is based on related calculations in [ DKM| . 



Lemma 4.2. Suppose dfi{z) = u(9)d9 is an absolutely continuous probability mea- 
sure on the unit circle and u(9) = u{—6). Then g(z) has the following properties 

1. g is analytic in C — S U (-co, —1). 

2. On (— oo, — 1), g+{z) — g-(z) = 2ni. 

3. g(z) = logz + 0(|) as z — > oo. 

4. e q9 ^ is analytic in C — S. 

5. e q9( - z ~> = z"(l + 0{\)) as z -> oo. 

6. (?(0) = ni. 

7. g+(z) + g-{z) = 2 J_ log |z — s\dfJ,(s) + i{4> + n) on z G E where <j) = arg(z). 

8. .g+(z) — g-(z) = 2ni f£ d/i(s) on 26 E. 

Proof, (i)-(v) are trivial. For (vi), 

pit pTT 

g{0) = / log (0 - e l6 )u(9)d9 = i{9 + n)u{9)d6 = ni 



using the evenness of u{9). 

For (vii), fix z = e 1 ^ £ E. Then arg(z — e ) is analytic if — n < 9 < <fi and 



<?+(z) = / log \z — e lS \dn{z) + i / arg(z —- e ie )d[i(z) + i / arg + (z — e ie )d/i(z), 

J —TT J —TT J (f> 

/TT /°<p I°TT 

log \z — e ld \d[i(z) + i / arg(z — e t0 )dfi(z) + i / arg-(z — e i9 )dfi(z). 
-TT J —TT J 4> 

Note that for tf> < 9 < tt, 

arg+ie** - e lB ) - arg-ie* - e 10 ) = 2tt. 

This yields 

/TT pTT pTT 

log\z-e ie \dfi(z)+2i / arg + {e l(t, -e ie )d^{z)-i / 2nd i i{z). 
-7T J —TV J <p 

Set 

/»7T /»7T 

F{<j>) = 2 arg(e l<l> - e l9 )dfi(z) + 2 / arg+ie* - e l9 )dn{z) - 2tt \ d\i{z) - cp 



If we show F(<f>) = tt, then (vii) is proved. Note that (a) ar<?(e l * — e % ^~ ) = <f)+ 
(b) arg+ie* - e^+) = 0+ ^ and (c) -harg(e^ - e ie ) = i. This gives us F'(<p) = 
2arg{e i ^-e l4 -)u{(j))~2arg + {e i ^-e l4 +)u{(j))+2TTu{(j)) = 0. But F(tt) = tt. Therefore 

F(<j>) = TT 
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For (viii), 

9+ {z) - g_{z) = i ["[arg+ie* - e w ) - arg-ie* - e i6 )]d^(z) 
J<f> 

r 

— if 2ndn(z). 



□ 



Let M be the set of probability measures on E. The equilibrium measure d/j>v(z) 
for potential V(z) — —%(z + z^ 1 ) on the unit circle is defined by the following 
minimization problem, 



(4.6) 



inf 



log|z — w\ 1 dfj,(z)dfj,(w) + / V(z)dfi(z). 



The infimum is achieved uniquely (see, e.g. |ST|) at the equilibrium measure. Let J 
denote the support of d\iy. The equilibrium measure and its support are uniquely 
determined by the following Euler-Lagrange variational conditions : 



(4.7) 



there exits a real constant I such that, 

2 / log \z - s\dfJ,v(s) - V(z) + 1 = for z e J, 

2 / log \z - s\dn v (s) - V(z) + l<0 for z e E - J. 

JT, 



In Lemma 4.3 below, we find d/iy, its support and I explicitly from this variational 
condition with the aid of Lemma 4.2. Let 



(4.8) 



g(z) = gv{z) = I log(z - s)dnv(s) 



where dfiy is the equilibrium measure. Following [DKMVZ1], we define 

to«(z) ee e^' 73 Y(z)e- q9 ^ 
Then nrS 1 ^ solves the following new RHP, 



(4.9) 



where = 



mS^{z) is analytic in C — E, 

(z) = m_ (z)v^ on E, 
mW(z)=J + 0(i) as z -> oo 



M9--9+) JL P q{g++g--v+l) 
e q(g+-g-) 



and 



(4.10) 



v g-l 



-Y 21 (0;q) = -rr^(0)e ql e 9a ^ = -(-l) g m ( ^(0)e ql , 



from Lemma 4.2 (vi) 



Once again we note that this RHP for is equivalent to the RHP for Y in 
the sense that a solution of one RHP yields a solution of the other RHP, and vice 
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versa. Using Lemma 4.2, the jump matrix i>W is given by 



(4.11) 



inside the support of d/iy. 



(-1)9 



\ er 

outside the support of d/j,y, 

/ e - 2 -3«/;<iMv(e) f^q^l J^g\z-e ie \dnv(e)-V(z)+l] 







As indicated in the Introduction, the purpose of the g-function is to turn exponen- 
tially growing terms in the jump matrix for the RHP, into oscillatory or exponen- 
tially decaying terms : this can be seen explicitly in (4.11), using (4.7). 

We have explicit formulae for the equilibrium measure and I. For < 7 < 1, the 
equilibrium measure has the whole circle as its support but for 7 > 1, a gap opens 



up. See, also [GW] and |Jol 



Notation : xb{&) denotes the indicator function of the set ficE. 



Lemma 4.3. For the weight V(z) 
are given as follows : 
1- IfO < 7 < 1, then 



;(z + Z ), the equilibrium measure and I 



(4.12) 



d/j, v (0) = — (1 +7 cos 8)d8 
2ir 



and 1 = 0. 
2. 7/7 > 1, then 



(4.13) 



dfiv(O) = -cos(-) 1 



- - ^ 2 {-)x[-e c: 9 c ](0)d0 



(4.14) 



where sin 



condition (4.7) is strict. 



I = — 7 + log7 + 1, 
-,0 < 9 C < 7T. In this case, the inequality in the variational 



Proof, (i) First, it is easy to check that djj,y{&) defined above in ( 4.12 ) is a positive 
probability measure. We set 

log(z - e l6 )^-(l + 7 cos 6)d6. 



9(z) 



2tt 



Then 



/ / , 1 f 1 7 / 1 n n ds 

9(z) = — / 1 + + -• 

Zm J s z — s I s 



Using a residue calculation with g(z) = logz + O(j) as z 
i W log e i{ - 0+v \l + 7 cos 9)d0 = iri, we have 

'logz-^ \z\ > l,z<£ (-oo,-l), 

— + 7TI \z\ < 1. 



00 and g(0) 



(4.15) ff (z) 
Therefore we have 



9+{z) +9-{z)= log z - -{z + z 1 )+TTi. 
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From Lemma |4.2| (vii) , we have 



2 / \og\z-e 10 \^-(l + jcos9)d9 + ^-(z + z 
J-tt 2tt 2 



for any z = e 1 ^ with I — as log z = i<f>. 

(ii) It is straightforward to check that the above measure ( 4.13Q is a positive 
probability measure. For g(z) defined in as before, we have 



g'(z) = [' -^2 C os(fU±-sin 2 f^ 
J_g c z - e %0 7r 2 y 7 2 

= T~ 1^ 1 S+ 2 V ( a -0(«-f)^ 

where £ = e 10 ", and the branch is chosen to be analytic in C — {e 10 : 9 C < \9\ < ir} 
and v(s — £)(s — f) > for real s >> 0. From a residue calculation, we obtain 

Integrating, we have for |z| > 1, z ^ (— oo, —1), 



= 1 iogz - 2(z + z- 1 ) + l + ^y/^-Ois-i)^ + 



and for |z| < 1, z <^ (-1, 0], 



s (*) = ±iogz- + + 2 + f ^ v /( s -0( S -0- + . 9+ (i) 

2 4 2 4 y 1 _ s v sz 

where <?+, g_ denote the limit from inside and outside each and l+o, l_o denote 
the outside and inside limits, 
(a) For H < C , 

9+{z) + g-{z) = log* - |(z + z- 1 ) + 7 + + fl-(l) 

From Lemma |4.2| (vii) , we obtain 



/ 9 I \ 9 
log 1 1 - e ie | ± cos - J - - sin 2 -d9 + i(0 + tt) 
-e c tt 2 y 7 2 

= — [ log(2| sin ^|) cos ^W- - sin 2 ^d6» + 7ri 
7T Jg 2 2 y 7 2 

= 2 log 2 - log 7 + - / Iog(sin0)cos 2 0<f0 + 7rz 
77 Jo 

4 /"^ 

= 21og2 - log7 - 1 + - / Iog(sin0)d0 + ni. 
f Jo 
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after simple change of variables and integration by parts. But we have J Q 2 \og(sin9)d9 
— -| log 2 from 

2 / 2 log(sin 9)d9 = I * Iog(sin0)d0 + / Iog(cos0)e20 
Jo Jo Jo 

Iog(isin20)e20 



= --log2 + /Iog(sin20)d0 + / " Iog(sin20)c20 
2 Jo -'f 

= -^log2 + / 2 Iog(sin0)d0. 
2 Jo 



Therefore 



9+i z ) + 9-( z ) = logz - -(z + z ) + 7 - log'7 + 



From Lemma |4.2| (vii), we obtain the desired result for \<j>\ < 9 C , z = e 
2 [ log \z - e ie \dn v (9) + l{z + z~ x ) -7 + log 7 +1 = 0. 

J-7T ^ 

(b) for 9 C < 4> < tt (— 7r < < — C case is similar), 

+ = logz - ^(z + z- 1 ) + + ,g_(l) 

7 + 1 



i(f> 



But 



c .2 



Therefore, using Lemma |4.2| (vii) and calculations in (a), we obtain for |0| > 9 C 
2 / log \z - e ie \d f i v {9) + ^(z + z- 1 ) - 7 + log 7 + 1< 0. 



□ 



In the following Sections, we distinguish the two cases, 7 < 1 and 7 > 1, due to 
the difference of the supports of their equilibrium measures. 



5. < 7 < 1 

From ( 4.15| ) , we have the explicit formula for the g- function 



9(z) 



logz-^ \z\ > l,z <£ (-oo,-l) 



-32! + iri \z\ < 1. 
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With this <?, I = from Lemma |L3J (i), and our RHP ( |4.9| ), or equivalently (4.11), 
becomes 



(5.1) 



(i) 



to' 1 ' analytic in C — E, 

/'(-l)«««e^(*- a ~ 1 5 (-l)? 

(-l)9z-«e^ a ^- s,_1 ) 
= J + 0(|) as z -> oo 

and k 2 ^ = _(_i)g m W(o) f rom ( p^| ). 
We define to' 2 ' in terms of to' 1 ' as follows : 
for even q, 

{to' 2 ' = to' 1 ' 
m (2)= TO (l)(0-l) 

for odd q, 

U< 2 > = (S-i)™ (1) U-i) |z|>1, 
\to' 2 ' S ( 1 _« 1 )to( 1 )(_ 1 - 1 ) N<1- 
Then we have a new equivalent RHP 



(5.2) 



|*| >1, 
\z\ < 1. 



,(2) 



(5.3) 

where w' 2 ' 



i (2 V 2 ) on E, 



4 ' 2 ' =/ + 0(i) as 



(-l)«z-«e" 



-(-l) q z q e 3 ^ ( - z ^ z 




on E, 



and k 2 j = 77122 (0). 




Figure 8. E' 3 ' and ft' 3 ' 



Introduce E' 3 ' = S( 3 '( 7 ) = E^ } U E^ t (see Figure |) as follows. For fixed 

7T/2 < |6>[ < 7T, 

(5.4) F(p) = F(p, 6) = Re(l{z - z" 1 ) + log z) = l{p- p' 1 ) cos 6 + logp, 
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where z = pe 
(5.5) 



9 < p < 1, has the minimum at 
1 



9 - 



7 2 cos 2 i 



-7 cos 9 



and F(p e ) < 0. ( Note that p e < for |0| < tt/2. ) For § < 7 < 1, we take 



(5.6) 



,(3) 
J in 

,(3) _ 



{pee 19 : 3^/4 < |0| < vr} U {p 3 „ /4 e ie : |0| < 3^/4}, 



= W^ 9 : 37T/4 < |0| < tt} U {p^e* : \0\ < 3^/4}. 

Orient £( 3 ) as in Figure |[ And finally, for < 7 < ±, set £ (3) (7) = £( 3 )(i). 
Of course, 

E (3) 

varies with 7 S [0,1]. However, using estimates from |GI 
is not difficult to show that the Cauchy operators C± on L 2 (£( 3 )) are uniformly 
bounded, 



(5.7) 



|C 4 



y ±||L2( S (3))_^ L 2( S (3)) 

for all < 7 < 1. Observe also that in the limit 7 
cross 

±\x + 11 for z = 



< C < 00 

1, £( 3 ) takes the form of the 



-1. 



(5.8) y — ±\x + l| tor z = x + ly near 

Apart from the neighborhood of z = — 1, there is considerable freedom in the 
choice of E^ 3 ^. For example, 37r/4 could be replaced by any angle between it/2 and 
n. Also the form of the contour for \8\ < 37r/4 is not critical, as long as it has the 
general shape drawn in Figure ^ : all that we really need is that the jump matrix 
u' 3 J below has the property supr zeS (3) : | ar g(z)|<37r/4} l w< ' 3 ' ) —^1^0 exponentially as 
q — ► 00. 

Using the factorization 



,( 2 ) 



we define 



1 







(-l^z-vHPt 2 - 2 ') 1 







r m (3) = m (2)( 6 p))-i 

m (3) = m (a)(6(_ 2 ))-i 

m (3) _ m (2) 

Then m^ 3 ) solves the RHP (£( 3 \u( 3 )) where 

'l -(-l)9z«e¥( z - z_1 ) 



in Q 
in £7 



in n< 8 \fi< 3) . 



(3) 

2 J 
(3) 

3 i 



(5.9) 

and 
(5.10) 



„(3) 



,(3) 




on S 



(3) 



on £ 



(3) 



<i = 4 3 2 } (o) 



As g -> 00, w< 3 )(z) -> I. Set £°° = £< 3 ). The RHP 
(5.11) 



m oo =m ooj Qn SOOj 

m°°=J + 0(i) as z^oo 
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has, of course, the unique solution, m°°{z) = I. 

Let < 7 < 1 — 5\ for some < <$i < 1. From the choice of E^ 3 \ 



ll« (3) -/|| 



(5.12) 



L°°(£( 3 )) = Sup 

37r/4<e<57r/4 



\e F{pe > d) \< sup I. 

37r/4<9<5ir/4 



But, for < 7 < 1, a straightforward estimate shows that 



(5.13) 

so that 
(5.14) 



2^2, 



yT _ ^ + log(l - VT - - log 7 < -^-(1 - 7 ) 3/2 



k (3) ||oo 



h) (3) -JL<e — ~ Al 9 



as q — ► oo. 



Since HC^cs) || L 2^ L 2 < C||w/ 3 ) ||oo, for some constant C independent of 7 (see (5.7)), 
(/ — C w (3) ) _1 is i nvertible for large q and the solution for the RHP (S^ 3 \ z/ 3 )) is 
given by (see (|2.9|)) 



(5.15) to (3) (z) = / + 



S(3) 



and (see ( |5~10| )) 

(5.16) ^_ i=m (3) (Q) = 1 + 



((/- cu )- 1 )/)( 5 )(^ 3 ) («)-/) d S 



2?ri 



Now from the fact that the length of T,^ is uniformly bounded and dist(0, S) > 
c > for all 7 £ [0, 1], we obtain, 



(5.17) 



\^ 1 - 1 \< C W vi3) - I \\oo<Ce-^ /2 i 



The above calculation also applies to the case when 7 — > 1 slowly. Indeed, 
suppose I < 7 < 1 — 3T7§ 273 ' wnere -Wi > is a fixed, sufficiently large number. 
(The lower bound h is chosen for convenience. Any fixed number between and 1 
would work.) From (5.13), (5.14), for some constant C which is independent of 7, 



(5.18) 



(3) ||L2^ L 2 



. 2-1/2 „Ci_„-|3/2 ^ 2V2 J,/ 3 / 2 , 1 



- <2<1, 



if Mi is sufficiently large. For convenience, we only consider Mi > 1. From (5.15), 



,(3) - 



1 



[(J-C^,,)- 1 ^*)^*) 



ds 



= I 



27T« ./j](3) s — z 

I /• u;( 3 )( S ) + [(/-C„ ( 3 ) )- 1 C,„ ( 3 ) /]( S ) W ( 3 )( S ) 



2™ 



S( 3 ) 



s — z 



and, as diag{w^ > ) = 0, 



(5.19) ^ = m| ) (0) = 1+ (J- ( [(/-C^ca,)- 1 ^,/]^)^ 3 )^^) 



SO) 
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(5.20) 



Hence, we have 

K_! - 1| < W-c^^c^iuq—^-u* 

H 2ms 

<\\(i-c w w)-^ L ^ L 4c wi sj\\ L 4—^\\ L z 

I-KIS 

<C\\w^\\ L ~\\w^\\ L , 
<Ce^^ 3/2 \\w^\\ Ll , 

where the (final) constant C is independent of 7, q and Mi (sufficiently large), 
provided that < 7 < 1 — 2 17^573 ■ 

Since the length of E^ 3 * 1 is bounded, we have 

|M 3) II u <Ce-^^-i?'\ which 



is the same estimation (5.17) as in the case 7 < 1 — b\. But for future calculations 
(see (7.2) below), we need a sharper result. We estimate ||i// 3 '||£i as follows : Focus 
on E^. For E^ t , similar computations apply. Only the 12-component of u/ 3 ) is 
non-zero. Set 9 — -^75 \ogq. 

|«9e*(*- a_1) ||dz| = (l) + (2) 



where (1) is an integration over \8\ < n — and (2) covers the remainder. Note 
from O) that \dz\ < Cd8. Substituting p e O) into F(p,6) O), we obtain on 



v (3) 



(5.21) l^e^C*-* -1 )! < e «(\/ 1 -')' 2 cos2 e+io g (i-^/i- 7 2 cos 2 8) _ log( _ 7 cos e) ) ^ 

Setting, 7 — ► — 7COS6' in ( |5.13|) , we obtain for z £ Sj„ , 

(5.22) |^e^^ _I )| < e -^(i+7cose)^ < e -c q (n-\e\f 



Hence, adjusting the constants C if necessary, we have 

C_ 



(1) < Ce~ c ^ < 



and 

/•<? , /-log? „ ,7/ r7 



Therefore, 

C 



(5-23) || W (3)|| L1 <__ 

and we obtain 



(5.24) l^-x - II < -^e-^ 1 -^ 2 . 



Let M2 > be a fixed number and consider 1 — 2 1 / 3 q 2 / 3 — 7 — 1- For this 
case, as q — > 00, 7 — > 1 and pe =7r — > 1. We need to devote special attention to the 
neighborhood of z — —1, where we will introduce a parametrix for the RHP, which 
is related to the special solution of the Painlevc II (PII) equation ( |l.4| ) given in 
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Section 0. For a discussion of parametrices in RHP's, see e.g. [DZ2|, [DKMVZlj. 
Set 7 = 1— 2l/3 * 2/3 ■ The region above corresponds to < t < M%. Let O be a 
small neighborhood of size e around z = — 1, where e > is a fixed number which 
is small enough so that first, 



(5.25) 

and second, 
(5.26) 



the map u defined below is a bijection from 0, 



the inequality ( 5.3 2| ) below is satisfied. 

(3) 



The goal is to solve the RHP for m> ' explicitly in this small region. 

Let u = \{z — z^ 1 ) in O. As noted above, we choose and fix e > sufficiently 




Figure 9. 

small (in fact, any number < e < 1 would do.) so that z — > u(z) is a bijection 
from O onto some open neighborhood of in the u-plane : under the bijection, 
EnO becomes a part of the imaginary axis. Set 

qV*u(z) _ g 1/8 l .-i, 

~~ i2 4/3 -, 2 W Z ^ } ' 

Note that with e fixed, there are constants c\ , C2 > such that 

(5.27) ck? 1 / 3 < \\{z)\ < c 2 q x l\ 

for all z G 9C Under the map z — > A(z), EflO now becomes part of the real axis 
and 

(5.28) A(£< 3 > H O) = {x + iy : = ^ Tj^/ , \x\ < 

1 _i_ if / 7 x 
1 + 9 V3 

where c is a fixed small number. As §373(1 — I 2 ) = |(l — 2 4 /4 2/3 )' ^ — * — 
we see that the contour A(£( 3 ) fl C) makes an angle < 7r/4 and uniformly bounded 
away from zero as q — * cxd, hence has the general shape of the contour in Figure ^ 
Section || within the ball A(O). We define 

(5.29) S p// - 2 n A(O) = A(E (3) n O) 

and extend S P77,2 smoothly outside A(O) in such a way that it is asymptotic to 
straight lines making angles between and 7r/3 with the real axis, ft is clear from 
the estimation in Section ^, and the preceding calculations, that for such a contour 
£ p// > 2 , the bound ( gjlj ) for the solution m PII ' 2 (z,t) of (£ p// > 2 , vf 11 ' 2 ) 



M 2 

(5.30) is uniform for 7, q satisfying the relation 1 , . < 7 < 1. 

2 ' a 1 



LONGEST INCREASING SUBSEQUENCE 



2!) 



Introduce the parametrix around z = — 1 as follows. Define 



m p (z) =m ru ' 2 (X(z),t) 



0-E( 3 ), 



)m p (z)=I in c -£( 3 ). 
As q — > oo, |A(z)| — > oo for z G <90, and we have for v p (z) = v^' II ' 2 (\(z)), 
' ' m p {z) is analytic in C - (E (3) U dO) 



(5.31) 



m p+ (z) = m p -(z)v p (z) on OflE' 3 ' 
m p+ 0) = m p -(z)I on O c n £ (3) 

^to p+ (z) =7+ mi A(z) + 0( xct) on dO as <? -> oo. 

The key fact is that v p is an approximation to i/ 3 -* with error of order ^573. 
We compare, for example, the 12-components of v*- 3 -* and v p on Y,^ 3 K We focus on 
O n sj- ■ Using the it variable, the 12-entries of and u p are 



exp 



(q[yu + log(\/ 1 + u 2 - it)]) 



and 



■ exp(g[7u — u + — u 3 ]) 



respectively. By ( |5.2lD and ( |5.22[ ), we have for z£0n s£\ 

| e g[7«+log(\/T+^-«)]| = | z ? e ¥( 2 - z_1 )| < e - ^ (1+7 cos S) 3 /^ 



From the Taylor expansion of the odd function log(Vl + u 2 — it), 

log(\/l + u 2 — it) = —u H — it 3 + u 5 r(u), 

6 



40 



where r(0) = 

sup| u |<i |r(tt)|. Note that for 



and r(it) is bounded for small u, say |u| < i. Set c = 
= p e e ie e On S^ 5 (see §J)), as g -> 00, we 



have |m| < c(l + 7COS0) 1 / 2 < c^'e) 1 / 2 . Therefore, if we have chosen e > small 
enough so that 



(5.32) 

we obtain, as q — > 00, 



2\/2 A „ 5 , ^ 1 



| e g[7«+log(Vl+« 2 -«)] _ e «[7«-«+H" 3 ] I 



q[7U+log(v / l+M 2 -M)] I 



-qu 5 r(u) I 



(5.33) 



< e-^^« e ) 3/2 x g | U | 5 ||r|| L . ({|u| <, })e ^) 5 (^- «) 5/2 

< C g (l + 7 cos0) 5 / 2 e ^-^ 2+ ^ 5c ' e )( 1 +^ osf, ) 3/2 
<C g (l + 7 cos0) 5 / 2 e-5(i+7co S ^/ 2 

C 



< 



,2/3- 
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where we have used the basic inequality |1— e z \ < |z|e' z ' and the fact that \\x 5 ^ 3 e~ x 
C. Since 

'\ — e lh u + lo S(Vl+u 2 -u)] _|_ e q[7ti-u+i« 3 ] 
,0 1 



3/2 



V&V- 1 = 



on ons 



(3) 



we have 



1 



(3) 

For O fl So Ut , we have a similar estimation. On the other hand, for O c , the error is 
exponentially small ; \\v^y- 1 - I\\ L ^ (0 -n^)) = ll« (3) ~ ^IL«>(e)<=ns(3)) = 0(e~ C9 ). 

Now define i?(z) = m^m" 1 . The ratio is analytic in C — (E^ 3 ' U <90) and the 
above calculations show that the jump matrix vr = mp-V^'Vp^m^ satisfies 



\VR - I\\oo < g 2/3 



(5.34) { IK-/IU < C*e- C « 
Vr = Up 1 = m"i = I 



on o c ns< 3 ', 

MT) — ' °M*(-jfas) on dO, as g -> oo. 

In (5.34), we have used the fact that m PII,2 (z,t), and hence m p , is invertible and 
bounded for (z,t) E C x [0,M 2 ]. (see, (^30|)). 

From ( H|4j ) and ( ggg ), we see that I||oo = ||iur||oo < C ^iff > ■ In particular, 

(7 — C^jj) is invertible for large q and by (2.9), R is given by 



R(z) =1 



2ni 



n(s)(v R - I) 



ds 



£< 3 >uae> s — z 



where fl solves (I — C WR )fi = I. As \\v R — 1]]^ < Ct ^J^ , we have — 7|| L 
Cm 2 (-^73), and also 



i? 22 (0) = 1 + 



2m 

C(M 2 



£< 3 >uae> 



ds 1 
(vr - l) 2 2(s)— +Om 2 (^tt)- 



Thus, using — Jjloo < 2/3 2j in (5.34) for the second equality, and m 1 



PII.2 



f 77 ' 1 (see ( [2.20 )) for the last equality we obtain 

K^-l = ^22(0) 

= l + ^-l (^-1)22(2)— + o M2 (-L) 

m i.22 0) dz 



1 

2?ri 



A(z) z 



+ °A^(-273") 



(5.35) 



m l,22 W 



1 dtt 1 

2m J u{90} —§± u (_ V^TT) + M2 W 



i2 4/3 pijo, , „ , 1 . 

^1,22' W+0 A / 2 (^) 



z2 4 / 



= 1 + TT7r TO f 2 2 4 W + %(^) 



1 



9' 



Note that error in (5.35) is uniform for < t < Mi- 
We summarize as follows. 
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Lemma 5.1. Let M\ > be a fixed number which is sufficiently large so that (5.18) 
is satisfied. Also let M% > and < S\ < 1 be fixed numbers. As q — > oo, we have 
the following results. 

1. // < 7 < 1 — S\, then, for some constants C, c which may depend on Si, 

K-i - 1| < Ce~^. 

2. If ^ < 7 < 1 — 2l ff \ /3 ; then, for some constant C which is independent of 
Mi satisfying ( 5.18 ), 

2 . _1|<- JL p -^#9(1-7) 3/2 



u 2 _i-^!rti<w 

1 g l/3 m l,22l E J| ^ g2 /3 • 

where t is defined by 7 = 1 — 21/3 ^ 2/3 . 



6. 7 > 1 



Let 9 C be given in Lemma |4.3|, sin 



, < 9 C < it. Decompose £ = Ci U C 2 



where Ci = {e J 



< |0| < tt} and C 2 = £ 



Ci. Note that on the support 
>.](0)d0 = 



of the measure d^iy in ( 4.13 ), d^v{0) — ^ cos(|)y i — sin 2 (|)x[ 
^^V^-O^-C- 1 )^ for * = e". 




Lemma 6.1. Define a(z) — — j ^r-\/(s — £)(s — £ _1 )ds where £ = e l6c and 

the branch is chosen to be analytic in C — Ci and ^/(s — £)(s — > /or real 
s > 0. TTien 

1. e 2a is independent of the path in C — ((7i U {0}),. 

2. exp(-27ri = cxp(2a(z)) /or z = e^, |^| < 6> c . 

3. exp(2/_ r 7r log|z-e ie |^y(6')-y(z) + /) = exp(-2a_(z)) for z = e^, |0| > 

0c- 
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Proof. Property (i) follows from a standard reside calculation : the change in a(z) 
around the point at is — ni, and the change in a(z) around G\ is 0. Property (ii) 
follows from the definition of a{z). For (iii), set 

F(4>) = 2 f log|e # -e 4e |^v(6»)+7C0S(/) + / + 2a_(e^) 



for z — e 1 ^, |0| > 6 C . From the variational condition (4.7), we have F(8 C 
Differentiating, 

2e^ 



L e^ — e' 
Z7 + 1 



l\dnv{0) - — ie - < 



iJ4 _ p -i<t>^ 



Y / (e^-C)(e^-C- 1 ) 



j_ 

2tt 



z s + 1 



V(8-0(3-£- l )d S 



c 2 



7 



V) z + 1 



^-otz-e- 1 ). 



2i v 2 z 

A residue calculation similar to that in (i), now shows that F'((f>) = 0. Therefore 
we have F{<j>) = 0. □ 



Note that e 



1 for <p outside the support of d[iy , i-e. for \<f>\ > 8 C 



By ( 4.11 ) and above Lemma, our RHP becomes 
'm^(z) is analytic in C — X, 



(6.1) 



,(1) e 



-2<ja 





(-1)' 



(i) (i) / 1 (-l) ? e- 2 ^- 
ml = m 
+ " ^0 1 

= J + O(-) as z -> oo 



on C 2 , 
on Ci, 



and k 2 ^ = -(-l)'?m^ 1 1 ) (0)e'?' = -(-l^e^+'^+^m^ (0) by fl4TC| ) and (p^) . 

We use the same conjugation ( |5.2[ ) for as in the case 7 < 1. Then our new 
jump matrices for rrS 2 ^ are 



(6.2) 



,(2) 



1 -e" 2 ^ 
, e 2, ? a 

e -2q Q _ A 

1 j 



on C2 
on Ci 



and n 2 q _ x = e'^-T+^T+^m^^O). 

Set E (3) = Ci U Cinside U C ou tside where C insid e and C outS ide are open arcs as 



chosen below. Note the factorization i/ 2 ) = 



1 

e 2?a X 



_ e -2ga 
1 



on C2. Set 



Rea = i?, /ma = I so that a — R + ii . Recall the Cauchy-Riemann equations in 
polar coordinates (r, 0), 



<9i? 9/ 97 OR 
dr d9 ' <9r 88 
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For z = e l9 S C%, a(z) = — 7ri f„ d[Av{8') is pure imaginary and 



dl 9 , ■ \ 7 A 1 ■ 2A n 

7^ = lisH") = ~^-cos(-)J - -sm (-) < 0. 



d9 86 



7 



Hence 



np or 

i? = and — = — <0 on C 2 . 
or o0 



Therefore for fixed 9, e %e £ C 2 , there is ei = ei(0) > such that R = Rea > 
(resp. < 0) for z = re 10 with 1 — e\ < r < 1 (resp. 1 < r < 1 + ei). We take Ci ns id e 
(resp, Coutside) such that |e~ 2a | < 1 (resp, \e 2a \ < 1) on Cinside (resp, C7 out sirfe)- 
Clearly there exist Q < pi,p2 < 1 such that |e~ 2Q | < pi (resp, |e 2a | < pa), for all 
z G Cinside (resp, z € Coutside), apart from a small neighborhood of the endpoints. 
Introduce the regions , k = 1, 2, 3, 4 as in Figure fill Define to' 3 ' as follows, 



m (3) 

m ( 3 ) 
m (3) 



TO (2) 

m (2) 
TO' 2 ' 



1 -e~ 2 « Q 
,0 1 

1 V 



in £7 



(3) 

3 i 



(3) 



Then v^ 3 ' is given by 



1 

.0 1 

1 

e 2qa 1 



in fi^fif. 



Oil Cinside -> 
Oil C outside j 

on Ci , 



and 
(6.3) 



v q-l 



- P 9(-7+log7 + l) rr) (3) 



m^(0). 



From Lemma 6.1 (iii) and the second variational condition in (4.7), we have, for 
any z G C±, 

e~ 2qa - ^ as q -> oo. 



(Recall that the inequality in (4/7) is strict from Lemma 13 (ii)). Also from the 
choice of C ms ide and C ou tstde, 

e~ 2qa -> 0, e 2 « Q ^ as q -> oo 



on C ms ide, C ou tside, respectively. Therefore w (3) — > w°° as g — > oo, 

On Ci n side U Coutside, 



(6.4) 



lj 

-1 

1 / 

The following result can be verified by direct calculation. 



on C\. 



34 



JINHO BAIK, PERCY DEIFT, AND KURT JOHANSSON 



Lemma 6.2. RHP (Ci,v°°) can be solved explicitly. 
(6 5) m»-n i (/3 + /3 " 1) ? (/? " r 

where (3(z) = ( ^T-i ) 1 ^ 4 ; is analytic in C — C\ such that (3 



Hence we expect from (p.3|), k 



9 (- 7 +io g7 +i)_^ j because m%{0) = 



Our goal now is to show that indeed — > m°° as q — > oo. As in Section^, we 
must control the behavior of the solution of the RHP for mA 3 ' near the endpoints, 



where the rate of exponential convergence u' 



3 , becomes smaller and smaller. 



Let #3, Mi > be fi xed numbers, let < 84 < 1 be a fixed, sufficiently small 
number sa tisfyin g ( 6.35) be low, and let A/3 > be a fixed, sufficiently large number 
satisfying ( 6.32 ) and ( 6.39| ) below. We consider 3 cases for 7 : 



1. 1 + 83 < 7 



2. 1 



+ 5r^<7<l + ' 



1 < 7 < 1 + 



Calculations similar to those that are needed for the asymptotics of the orthog- 
onal polynomial on the real line (see, DKMVZ1|), show that for 7 > 1 + 63, 



(6.6) 



o<?(-7+log7+l) 



o(i)). 



The error is uniform for 1 + £3 < 7 < L for any fixed L < 00. However, we will not 
use this result, utilizing instead (stronger) estimates from [ Jol (see next section). 




7 



Figure 11. 

We consider case (hi). Set 7 = 1 + 2l/3 ' i/3 with < t < M4 and u = sinO c = 
y/'y — 1. In defining E^ 3 ) above, there is some freedom in the choice of Cinside 



and C ou tside- We make the following choice (see (3.11) below). Set x 



2V 3 



t(l + 21/3* 2/3 ) 2 ~ — t < as 1 ~~ * °°; an d let T, PI1 ' 3 be the contour defined i 



m 



Figure § for this specific x. Let E' = {u = ^\ : A e E p// < 3 } = Uf =1 E' fe! and let 



€ > 



be small and fixed (see (6.7), ( |6.10| ) below). For definiteness, we can, and do 



assume that the rays E^, . . . , E4 make an angle of 7r/6 with the real axis. Consider 
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O' = {u : \u\ < e'}. If q is la rge e nough, then uo G O' . Set u — u(z) — ^(z — z^ 1 ). 
We choose e' such that (cf. (5.25)) 

(6.7) u is a bijection from an open neighborhood of z = — 1 onto O' . 

Clearly there are constants Ci, Ci > such that ci < < c-i for all u G 9C. 

Under tt , the points Uq, — uo are mapped into £, £ respectively, and w _1 (£g) = 
C\. Consider a point z G u _1 (£ 4 n O'), the inverse image of a point u G £ 4 fl O, 
Changing variables twice, v — ^(s — s _1 ) and w = v 2 , 

(6.8) 

-2a(*) = £ f^+^V^ + ^^-^ + e- 1 )- 



-17 



1/2 

v(l + h(v 2 ))dv 



(tu - u 2 ) 1 ' 2 [1 + - h{ul)) + h{ul)} dw 



17 

= _ T 

_n f 2 , _ 2)1/2 [( fc ( u §) - - 2k(u 2 )] (1 + feH)rf w 
4 7„ e lw u ° j (1 + k(u 2 ) + fcH) 1/2 [1 + (1 + Hu 2 ) + fcH) 1/2 ] 

= ^(« 2 - ulf/ 2 + o(\u 2 - u 2 \^ + o(\u 2 - u 2 \^ 2 ui 

where \Ju 2 — Uq is defined to be analytic in C — [— Uo, Uo] and positive for real 
u > uq ; h{w) — ~ 1) which is analytic in \w\ < e' and h(0) = ; 

k(w) = — w — 1), which is also analytic in \w\ < e' and k(Q) = 0. Since £ 4 is 

a straight ray of angle at uq, Re(—i(u 2 — u^) 3 / 2 ) < —\\u 2 — «o| 3 / 2 , which yields 

(6.9) I exp(-2a(z))| < exp(-^|u 2 - u 2 | 3 / 2 ) < 1, 

provided e' is sufficiently small so that 
(6.10) 

0(\u 2 - u 2 Q \^ 2 ) + 0(\u 2 - u 2 \V 2 u 2 ) <c\u 2 - ul\ 3 ' 2 {\u 2 - u 2 Q \ + ul) < ^\u 2 - « 2 | 3 / 2 



for u G £4 fl O' , where the terms on the LHS are given in (3.8). The same choice 
of e' gives rise the same result for z G u~ 1 (Eg fl C), and also | exp(2a(z))| < 
exp(-^|u 2 -wg| 3 / 2 ) < 1 for z G u _1 ((E' 1 U £' 2 ) fl O'). We thus fix £( 3 ) by choosing 

(6.11) S( 3 )=ti- 1 (S'nO') inside O' , 

and extending it to a contour of the general shape C\ U Ci na id e U C ou t S ide as in 
Figure [To| . 
Define 

f,n p (z) = m™(^ u (z),x) in O-SG), 
L p (z) = / in O c -E( 3 ), 



(6.12) 
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where m PII ' 3 (z, x) solves the RHP of Painleve II equation given by (|J|) and ( gTgg ) 



Then m p solves the RHP on T,^ U dO in which the jump matrix v p (z) is given by 

(6.13) U ~ ze^nO c , 

(m p +(z) z e dO, 



where v PI1,3 is given in ( p. 26 ) 



We compare and v p . First, let z G £ (3) n O such that u(z) 6 Ejfl O' . 
The 12-entries of and v p are — exp(— 2qa(z)) and — exp(— 2ig PI1 (^jj u(z))) = 
— exp(— ^r(u 2 — Mq) 3 ^ 2 ), respectively. With e' chosen small as above, using 

Re(-i(u 2 -ul)V 2 )<-\\u 2 -ut\ 3 / 2 , 
2 4M 4 

W 2 < 



7 2/3 



||x 5 / 2 e- 3/2 || LO o :[0jOo) <C, 
7 — 1 < 



2 1 /3 g 2/3' 



we obtain from (S.8) 



1/3 

| e -2ga( Z ) _ e -2 l3 p "(2 37 ^«( 2 )) | 



(6.14) 



< e i?. e (-i2(n 2 -^) 3 / 2 )| e -29a(,)+ 2l ^"(^ M (,)) _ 

< Ce" Al« 2 -^| 3/2 [,| u 2 _ U 2 |3/2 ( | U 2 _ ,2| + u 2 + (7 _ 

C(M 4 ) 



< 



q 2/3 



In a similar manner, for z such that u(z) £ E' 3 fl C, same result holds and for 
z G Coutside H 0, the difference of the 21-entries of w'- 3 -' and w p satisfies |e 29 "( 2 ) — 

1/3 

e 2i 9 (^3«(*))| < por z e Ee( _ i(u 2 _ u 2)3/2j = _| M 2 _ u 2|3/2 

Again by (3.8), the difference of the 11-entries of and v p satisfies \^ 2 i a -( z ) — 



n PH(3_ 



(f473 "( z ))| < £iMA, Therefore, we have 



q 2/3 

(6.15) H^V 1 - IWl^uo) < ^ • 

Secondly, for z G S< 3 ) n O c , |z - £|, | z ~ £ | > c > implies exponential decay 
for e - 2qa ^ and e 2qa ^ for z € C* in;jJ d e n O c and z G Gourde n O c , respectively. 
Therefore we have 

(6.16) Wv^v- 1 - /IL-^no.) < Ce" *. 
Finally, for z G dO, \u(z)\ = e' and 



1/3 PH,3( \ i 

(6.17) m p+ (z) = m ™(| i7lU (z), = J + -U + M4 (^), 

24/3 U{Z) Q 

by ( |2~25| ). Here the error is uniformly for < x < 2 A/4. 
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Now as in the case 7 < 1, define R(z) = m^ 3 'm p 1 . Then the jump matrix for R 
is given by Vr = mp^v^v^m" I. From ( |6l5| ), (|t|) and (|j7|), and also ( gjg ), 
we have 



II rll s' t ->( m 4 

WWII — J llL°°(£(3)nO) S: -^273 



C(M 4 ) 



(6.18) 



\\ V R - ^llL=°(s< 3 >no c ) ^ C e cq 

v R = m;l = i- 2 /3 q T/k u( ' z) lx) + Mi {qm) on DO 



onOntf', 

on o c ns( 3 ' 



(? 1 /3 M ( Z ) 

As in ( 5.35 ) for the case 7 < 1, using rnf 22 = m i22' 3 + (ix 2 /8) from ( 2.27 ), 



m| ) (0)=i? 22 (0) 



(6.19) 



124/3 w^/ ( 1 \ 



1 + 



*2 4 / 3 



"ifg(a?)-^-(l 



rf 2 



2 l/3 g 2/3, 



q 



1/3 



i 2 1 

< 22 (Z) + 25/3g l/3 + 



Therefore, from 



|) using x = -t(l + 2l/ i q2/3 ) 2 = -t + M A^m) and tne fact 
that -imf '^(t) is bounded for < t < M4, (this follows, for example from ( [2.18 ) 
and the boundedness of u(x) = 2imf I 1 2 ; alternatively statements like (2.2?) are 
true also for all the x-derivatives of m PII ' 3 (z; x), etc.), 



(6.20) 

t 2 ^ _ „9(-7+log7+l) rrl ( 3 ) 



/.•" .,=(''' ' ~' m 22 (0) 



eg (-7+logT+l)(l + ! 2 ^! m WI( a; ) 



<1 



1/3 



t 2 „ , 1 , 

+ Um ^7273^ 



25/3^1/3 



1 - 



2 5/3 g l/3 



z2 4 / 3 



+ OmA-) 1(1 + ^Ta-mfg^) + -75-^ + Om 4 ("^) 



' 2 / -<22(^)+Om 1 (^) 



J/3 



j2 4 / 3 1 

— m ^(_*) + M4 (_). 



Finally we consider the case (ii), 1 



Ma 
2 1 /3 ? 2/3 

with jump matrix v^ 2 ' given by (3.2), as follows 



< 7 < 1 + ^4. We conjugate 



(6.21) 

Define 5(z) 



■ft 



-,(4) 



,(2) 



m< 4 ) = m( 2 )(_° 1 i) 



«| > 1, 

z\ < 1. 



7_ . 

J 47TI 



— -\/( s — 0( s — £ 1 )^ s i where a(z) is the same as a(z) 



in Lemma OA , but now we choose the branch so that y/Js — £)(s — £ is analytic 
in C — C2, and \/( s ~ £)( s ~ ~ +s as s — > 00. Then the jump matrix for 
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:s<> 



As in the case of a, we have \e < 1 for z g Ci and |e"( z ) 
Cmside U C 'outside- Therefore taking g — * oo, we have 



(6.23) 



This RHP can be solved explicitly as in Lemma p.2j, and we find 




(6.24) 

where (3(z) 



,(6,00) _ 



^(/3 + r 1 ) htf-t 1 ) 



< 1 for z g 



( ^cli ) 1 ^ 4 is now analytic in C — C2 and /3 ~ +1 as z 



>(0) - and k. 



o<?( — C+log7+l) J_ 
V7 



as g 



DO . 

00. 



From ( 6.24 ), we have 

Again, we need to construct parametrices around £ and £ _1 in order to prove that 
indeed — > ttt/ 5 ' 00 ). Note that detm^ 5 ' 00 ) = 1 : this follows either by direct 
calculation or by a general argument as detz/ 5 ' 00 ) = 1. 

Let E" be the contour E" = 1 U R + e 2n% / 3 U M+e 4 ™/ 3 shown in Figure O. Let 




Figure 13. E" and 



w _ e 27r l /3 and get c see | D z2 



(6.25) < 



Ai(s) 


Ai(tjj 2 s) 






Ai'(s) 


lu 2 Ai' (lu 2 s) . 


e~ 




Ai(s) 


Ai{ui 2 s) ' 


e~ 




Ai'(s) 


lo 2 Ai'{uj 2 s) , 






Ai(s) 


—u> 2 Ai(u)s)\ 


e~ 




Ai'(s) 


-Ai'(ws) J 






Ai(s) 


—ui 2 Ai(u)s)y 






Ai'(s) 


u>Ai'(ujs) 1 


e~ 


'a <T 3 , 



< args < , 
— < args < it, 



ir < args < ^p, 



4p < ar#s < 27r, 
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where Ai(s) is the Airy function. Then W satisfies the jump conditions 

'1 l^ 



(6.26) 



vfi + = i$r_ 



o i / 

1> 

1 o\ 
i i/ 



z g 



z e 



z G 



2wi / 3 ,M+e 4 "/ 3 . 



Let Of and Of be neighborhoods around £ and £ of s ize e" V 7 — 1 , respectively, 
where e" > is a small, fixed number chosen to satisfy ( |6.27 ), ( 6.29| ) below. Since 



r|£ — £| = - V7 — 1 > e'V7 — 1; an( i have no intersection, provided 



(6.27) 



< e" < 1. 



For definiteness we assume that <90f and <90f are oriented counterclockwise. In 
Of, a simple substitution shows that a(z) — |(z — £) 3 / 2 G(z), where G is analytic 
and = (7 - l) 3 / 4 e -*(f e <=+f ^ Here (z - £) 3/2 = |« - f |3/ 2 e SW«-0 and 

6> c - tt/2 < ar 5 (z - < 6 C + 3tt/2. Define 

(6.28) \(z)^(z~0(G(z)) 2 / 3 , 

where (G(z)) 2 / 3 is analytic in Of and (G(z)) 2 / 3 ->■ (7 - l)l/2 e -i(0cW2) as z £ 
Of course, A 3 / 2 = la. It is a simple calculus question to verify that we may choose 



£ (Y-l) 




•OH) 











/ 





Figure 14. 

e" sufficiently small so that 

z — > A(z) is a bijection from Of onto an open neighborhood 
of in the A-plane, of radius ~ (7 — 1). 



(6.29) 



Define S (5) n Of = {z G Of : A(z) G E"}. As in the construction in ]DZ2| , set 
(cf. (4.34) in pZ2| ) 



(6.30) 



1 -1 
-1 -1 



5Fe^|(«-f)(<?W) 



2/3 



and for z G Of — S^ 5 \ define the parametrix for m^ 5 ) by 
(6.31) TOp(z) = ^(z)*(<7 2/3 A(z))e 95(z)<T3 . 
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Then m p satisfies the same jump conditions on E^ 5 ) nOj as 



i m 



(5) 



m v -V 



(5) 



And if q becomes large, then for z g <90 ? , |q 2/3 A(z)| > cg 2/3 (7 - 1) > cM 3 /2 1/3 
Therefore, 



(6.32) 



if M3 is sufficiently large so that the leading terms dominate 



in the asymptotics for the Airy functions in ^j(q 2 ^ 3 A(z)) (see e.g.[AS|), 
then by the explicit choice of E(z) in ( 6.30| ), we find for z S dO^, 

(6.33) rn p+ = m ^H Z )(l + 0(^^)^ = m< B ->(z) (f + _\ )3/2 )) ■ 

Noting the symmetry m/ 5 ) = m( 5 )(z), define E^ 5 ^ n 0^ = E( 5 ) n0£, and for 
z£O r E^ 5 ) set m p (z) = m p (z). We now extend E^ n (C ? U Of) to E^ to have 
the same general shape as in Figure Finally, for z € C— (0£ U0j U E^ 5 '), define 
m p (z) = m( 5 '°°\ 

Set O = C>£ U 0|. Then # = m^m p 1 solves a RHP on E^ 5 ) U dO with the jump 

on (E( 5 )nO)UC 2 , 

1 e^ 2 ^ 



matrix Vs = m„„i/ 5 )?;„ "m 



(6.34) 




,(5,oo)n 



(5,00) 



1 on CinO c , 

m^nte U Coutside) ^ ^ j 



on (C„ 

Let ei > be a fixed, small number : for example, we may take e± = e' satisfy- 



ing ( |6.7| ), (6.10) above. Choose 64 sufficiently small so that 
(6.35) £,£e {z: \z + l\ <ei} 



for 1 < 7 < 1 + 54. By calculations similar to ( p.8[ ) and ( |6.9| ), i?e(2o!(z)) < 
-c(7 - 1) 3/2 for z € (Cinside U Coutside) H {|z + 1| < ei} fl £> c (in fact the estimate 
is true on the full set (C msi de U C outS4( j e ) n {|z + 1| < ei}) and also |e 2<?5 ^| < e~ cq 
for zE (C insi de U C ouiside ) n {\z + 1| > ex}. Thus, | e 2 « 5 ( 2 )| < e -<=?(7-l) 3/2 for 
2 € (Ci ns ide U C ou tside) H C . Also, by calculations similar to ( |6.8D and ( |6~9| ) again, 
i?e(-2<5(z)) < -0(7 - I) 3 / 2 for 2 e Ci H O c . Therefore we have L°° estimation 



(6.36) 



Furthermore, from calculations similar to ( 5.23| ), on Ci ns ide H C n (Jm(z) > 0}, 
using |it + u | > |w — uq| °n the integration contour for the second inequality, 



-2qa(z) 



\\dz\< 



(6.37) 



< 



{u=UQ-\-xe~~ z7v / 3 :x>c^'y — 1} 



< 



c-^/7 — 1 

c 

9(7-1)' 
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The same calculations apply to the other part of Ci ns id e O c and also to C ou t S ide H 



On the other hand, length(dO) < 



O c so tha t L^(c inHde uc ou t Bide )no°) - g( 7 -i) 

CV7 — 1 an d length(C\ n O c ) < CV7 — 1, and hence, by the above L°° estimates, 
\K ~ AWpo) < C/{q{l - 1)) and \\v R - /|Ui (Cin0 =) < - 1)). Thus 



(6.38) 



C 



\\vr- I\\v-$iV»uao) < g ( 7 _ i) 



Using the choice of 0^ and Of, direct calculation shows that m^ 5 '°°\ hence 
(nj( 5 '°°)) _1 ^ ag d e i m ( 5 ,°°) — are uniformly bounded for 7 in the region 1 + 
2l/ ^ 2/3 < 7 < 1 + <5 4 , for z S O c - S( 5 ). On the other hand, even though the 

contour E^UdO varies with 7 and g, the length of T,^UdO is uniformly bounded 
f° r 1 + 2 i/t^ 3 2 /3 < 7 < 1 + <54- Also a simple scaling argument shows that the Cauchy 
operators C± on L 2 (E^ 5 ^ U 90) are uniformly bounded for 1 + 2l ^ 2 / 3 < 7 < 1 + 64. 



Therefore, 



l|C™ s llL 2 (s< 5 )uao)- > L 2 (E( 5 )uao) ^ ^ll w /}llL~(£( 5 )uao) 

C 



< 



(6.39) 



9(7 - I) 3 / 2 



+ C e -c a (7-i) 



3/2 



< 



C 



M. 



2/3 



Ce 



-CM: 



3/2 



provided that M3 is sufficiently large. 
From (2.£) and (6.38), we have 



1^22(0) -1| = 



1 
2~7ri 



w~ R -[{I-G w -)- 1 G Wk I]{z)w k {z) 



ds 



/E(5)uao \ • / 22 

< C(\\w R \\ L l^(5) ud0 - ) + ||W|j||^2( S (5) uaC |-|) 

< c \\wR\\Li(s(VudO) > as H w flJlL-(E(=)uao) is bounded, 
C 



< 



9(7-1)' 



and 



l-Kax(0)| < 



C 



9(7-1)' 

Therefore, from m$(0) = -R 2 2(0)m^ 1 ' oo) (0) + E 2 i(Q)mf 1 ,oo) (0), we obtain 



(6.40) k»_ x 



<?(-7+log7+l 



o9(~7+log7+l) . 



V7 V 



9(7-1)' 



Note that this is consistent with the result (3.t) for case (i) where 7 — 1 > 83. 
Summarizing, we have proven the following results. 

Lemma 6.3. Let 83, M4 > be fixed numbers. Let 84 > be a fixed sufficiently 
small number satisfying (6.35), and M3 > be a fixed, sufficiently large number 
satisfying (3.32) and d6.39| ). As q — > 00, we have the following asymptotics. 
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1. Ifl + 6 3 < 7 , 



oli- 



^o^+i)J_( 1 + iK\ 
\fi V i J 



where the error is uniform for 1 + 64 < 7 < L for any fixed L < 00. 



2. If I 



Ms 

2 l/3 g 2/3 



< 7 < 1 + h 



q(-f+l 



/TV 

where the error is uniform in the region. 

3. Ifi<i<i + ^m, 

z2 4 / 3 



l K 9-l 1 „!/ 



C(M 4 ) 



9 



where t is defined by 7 = 1 



2 l/3 ? 2/3 > 



< t < Af 4 



Note that, comparing Lemma 6.3 (iii) with Lemma 5.1 (iii), we have same result 



everywhere in the region 1 



2 1 /3<j2/3 



< 7 < 1 



M 



(6.41) 



\ K q-l 1 



i2 4 / 3 PII , C(M) 

-173-^1,221*) I < 2/3 ' 



where i is defined by 7 = 1 — 2l/3 * 2 /3 1 an d M is any fixed positive number. 



Also note from Lemma |6.3j (ii), that as q — > 00, 



(6.42) 



logK,_i < <7(-7 + log7 + 1) + ^Z— fi 



where is independent of M3 and is fixed once $4 satisfying (5.35) is determined. 



7. Asymptotics of n (A) as n — > 00 



In this Section, using Lemmas 5.1 and p.3| , we obtain the large n behavior of 
</> n (A). In the following £5, 5q, 67 are fixed numbers between and 1, and M 5 , M 6 , M7 
are fixed and positive. These numbers are free apart from the following requirements 



(a) S e satisfies ( |6.35 ), 

(b) M 5 > 1 satisfies flsTi 



(c) -Me > 1 satisfies (5.18), and 



(d) M7 > 1 satisfies ( |6.32 ), ( 6.39| ) and condition (7.8) below. 

We consider the following five cases for A > and n : 

1. 0< 

2 i < 
A - 2 — 

3. 1- 

4. 1 + 



*g <i_ 

n+l — 
2V 7 ! <^ 1 
n+l — x 






2 1 /3(„ + l)2/3 




M e 


< 2^5 < 1 H 

— n+l — 




2!/3(„ + l)2/3 


2i/3( n+ i)2/3 


M 7 


< 2^X < 1 H 
— n+l — 


- 8 6 


2l/3(„+l)2/3 
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Consider case (i). For any k > n, < 1 — <5s- From Lemma 5.1 (i), we have 
as n — > oo, 

OO OO 

(7.1) |log 0„(A)| = |^log4(A)| < J2 Ce' ck < Ce~ cn . 



k—n 



k—n 



Consider case (ii). We split the sum into two pieces. 

oo 

log^ n (A) = ^log4(A) 

k—n 

= ^ log 4 (A) +^ log 4(A) 

(1) (2) 

where (1) and (2) represent the regions 

(1) n+ 1 < k+ 1 < iVX, 

(2) 4VA < k + 1. 

M 5 



For (1), 2 ' fe^ 1 _>' :/, i ]- ; 



From Lemma 5A (ii), for some constant C, 



independent of M5 satisfying (5.18), 



|logK 2 fe (A)| <C- 



.2v5(fc + l)(l-|^) 3/2 



(fc + 1) 



1/3 



Using the fact that f(x) = 
second inequality below, we have, as n — > 00, 



, -"l 1 2 ^ x ) / is monotone decreasing in the 



(7.2) 

|£log4(A)|<C£ 



-2#(fc+l)(l-l^?) 3/2 



(1) 



(1) 



(fc+1) 



1/3 



<C e 

Jn+l 

< C(2VA) 2 / 3 

< C(2VA) 2 / 3 

< c 



2V ^ x(l 2%/A -)3/2 (ix 

^T73 

4v^Ay 3 / 2 



c- 



.2^2( n+ l)(l_!^)3/2 



1/3 



g 3(1 + B )V2 



(" + !) 

1 



(n + 1) 

+ Ce -^ n+1 )( 1 -^) 3/2 



("+!) 



(l + y)V3 



2VT 



ds 

e — + Ce 2 



<Cexp(-±(„ + l)(l 



i( n +i)(i-f^) 3 / 2 



2v^I \3/2 



n + 1' 



We use the change of variable y — — 1 for the integral in the third line. The 
fifth inequality is obtained from the substitution s = \f\y 3 / 2 , and at the end, we 



have used ^ < 1. 
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For (2), jj^j < |. Therefore, from Lemma |5.1| (i), we have 

oo 

|53 log K* (A) I < E Ce- ck <Ce~ cn . 

(2) fc+l=[4\/A] 

Summing up the above two calculations, we have, for case (ii), 
(7.3) | log 0„ (A) | <Cexp^ c (n + l)(l-^-) 3/2 ^, 

as n — > oo. Note again that the constants C,c can be taken independent of M> 



5- 



Consider case (hi). Set 
(7.4) M-I- ' 



n+1 2V3( n + 1) 2 /3 

so that — M 6 < t < Mg. We divide the sum into three pieces, 

oo 

log<£ n (A) = E lo § K fc(A) 

= E log 4 (x) + Y, l °s4W + J2 lo s K ' ( A ) 

(1) (2) (3) 

where (1), (2) and (3) indicate the following regions : 

(1) n+l<fc + l<(n + l) + ( M *~ ^ (n + 1) 1/3 

(2) (» + !) + ^|t^(- + !) 1/3 < fc + 1 < |(» + 1) - 4f> + 1)V3 

(3) |(„ + l)-^(„+ 1)V3 < fc+1 . 



For (1), as n — > oo, 



6Af 6 2\/A 2M 6 

1 „,„ < — - — < 1 



2 1 /3(fc + l) 2/3 fc + 1 2 1 /3(fc + l) 2/3 ' 

Hence from ( |6.41 ) , we have as k > n — > oo , 



l°g4(A) = _^ f mf^(2i(fc + l)«(l- |^)) +0^^). 



(fc + 1) 
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This leads to 



X lo s4(A) 

(i) 

[(„+!)+ <^ 7 _l) (n+1) i/3] 



E 

fc+l=n+l 



- m «J(2*(*+l)*(l-£^l))+0 



(fc + l) 1/3 ' * + l J,J6V P/3- 

<^(2 1/ V/ 3 (l-^))dx + M6 ( ' ) 



(n+l)+ ^7,3" (n+l) 1/J - 9 4/3 9,/T 

T l/3 m l,22^ a; U , ~iu 6 V 1/3 

(n+1) X X " 

(M 6 -t) 

21/3 i2V3 m ^ r 2 l/3 (B + 1)2 /3 (1 + 2 * + . . . ) 

' V 3(n+l) 2/3 
(n+1) - ^(n + iy/W ds , n , 1 



1 _ ; — 7T-, . ^1/3 ) m 1 » - — v + ° M ° (: 



(Me-') 

2 V3 



(n + l) + S (n + l) 1/3 V (1 + ^273 +■••) ' 6 ^ 1/3; 

«<MS ((« + - ^ + •■■>>- ^ + • ■ ■ >* 



+ ° M «(^I73") 



(Mg-t 
2 l/3 



2im[ I 2 I 2 (y)dy + Me 



1/3, 



,1/3, 



The fourth equation is obtained using the change of variable x = (n + 1) + 
s(n + 1) , and for the sixth equation, we use the fact that-^mf 22(f) is uni- 
formly bounded for —Me < t < Mq (see the remark below ( 3.19[ )). To pass from 
the second to the third line, note that for integers b > a, 



(7.5) 



6-1 .b b-l 

5^ /(»)-/ f(x)\<^2 s M\f^)-f(P)\--n<a,P<n+l} 

nr. n >J CL „ ,, 



< \\f'\\L°°( a , b )(b-a). 
For the case at hand, a simple calculation shows that 

Also, the contribution to the integral from the interval ((n+l)+ ^1/^ (n + l) 1 ^ 3 , [("■+ 
l) + ^f7^(^ + l) 1/3 ] + l) is M6 (l/n 1 / 3 ). 

For (2), 

1 2^ ±M 6 

- < < 1 TT7TT TT7T k > 71 — > 00. 

2 - jfe + 1 ~ 2V3(fc+ 1)2/3 
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As n — > oo, by a calculation similar to the case (ii), again using the monotonicity 
of f(x) = —^j S e~ 1 ^' x< - 1 ~ 1: ^^ 12 for the second inequality, we have 



|5>g4(A)|<C7£ 

(2) 



-2#(fc+i)(i-|r?) 3/2 

H73 



(2) + 

(n+i)-^(n+i) I/3 



(M B -t) fn , ,,1/3 



✓A-,3/2 dx 



< C(2\/A) 



2/3 



(n+1) 



WgX / H J U/2 



(. 1+5/ -> 



a.i/3 
r/i/ 



Ce 1 2 ; 



<C(2VA) 2 / 3 i/3 e 



<C/ _ds + Ce-^ 



(l+y)V3 
-^ 3/2 dy + Ce-^) 3/2 

"6 N3/2 



_/Wes3/2 

+ Ce 1 2 1 



1 i.r3/2 -iw g \ 3 / 2 



Mg. \3/2 



- A/ 



The first inequality follows from Lemma 5.1 (ii) (note that, by assumption, 
satisfies (5.18)). For the second line, in order to control the contribution to the 
integral from the interval [[(n+1) + ^0±(n + 1) 1/3 ], (n+1) + ^jfi(n + 1) 1/3 ], 

we use the inequality 1 + ^1757^77^273 

< 1 



^+1)2% < t for large enough n. For 
1, and for the fourth line, 



2V3( n+ l)2/3 

the third line, we use the change of variable y = 



we use the inequality prj^g > 1 — 2 1 / 3 (n+i) 2 / 3 — H ^ or sumc i cn tly large n. The 
fifth equation is obtained from the substitution s = \/\y 3 / 2 , and for the sixth line, 



we have used the inequality |^ < 1 

2VA 



Me 

2 l/3( n+ l)2/3 



< 2 for sufficiently large n. 



For (3), as tl — > 00 , -j^j < |, which yields, from Lemma |5.l| (i), 



1^ log 4(A) I <Ce~ cn . 

(3) 

Summing up all these calculations, for = 1 ; , * — 7577 

o r > „_|_i 2 1 / 3 (n+l) 2/J 

Mq, we have, as n — ► 00, 



with -M fi < t < 



I log 0™(A) 



2z<^ (2/ )^|<^§l + Ce-i^ /2 , 



for a constant C(Mg) which depends on Me, and for a constant C which is inde- 
pendent of Mq. Using the asymptotics of (see j2~T^)), we 
have (recall ~M 6 > 1) 



(7.6) 



I log (j) n {\) 



^ m i.22{y)dy\ < 



g(M 6 ; 



1 »^ 3 / 2 
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M, 



Now we consider case (iv), 1 + 2 i /3(ri+1) 2/3 - 

log n ( A) = ^ !og 4 + lo S K l > 

(i) (2) 



< < 1 + <V We write 



where (1), (2) indicate the following regions 



(1) n+l<k + l<2V\~ ^(n + l) 1 / 3 

(2) 2v / A-|^(n + l) 1 /3< A; + i. 



m 7 



For (1) , we have for n sufficiently large, 1+ 2i/3(fc+i) 2 / 3 — fc+i 
using ( |6.42 ), we obtain 

(7.7) 

E lo g K '( A ) 



< 2\/a < 1+ x Therefore, 



(i) 

fc+l=n+l v 

< — / xl f J ax 

! Ml 

• -i / 2l/3< " + " " !2x Ar '" " dy + 



C* 



2VA - (fc + 1) 

2^-^T(«+l) 1/3 C # 
n+1 2\/A — 

(2 VA) 2 + C* log f 2 f-;;+ 1} 

+i y \ (n + I) 1 /- 1 



-dx + C(M 7 ) 



C(M 7 ) 



i 

< -- 
~ 4 



«7 n + 1 

2 1 /3(„+i)2/3 2VA 



+ 1 
2VT 



(2VA) 2 (l- 2/ ) 2 d 2/ + C# 



2VA- (n+1) 
(n + l)V3 



C(M 7 ) 



< 1 [^^^ (2>/X) + C*(n + 1)^(^ - 1) + C (Mr) 
4j 1 _«+i n + 1 

2VA 



< -I( 2 \/a) 2 
- i2 V ; 



M 7 



2V3( n + 1)2/3 



2VA 



< 



+ C # (n + l) 2 / 3 (^--l) + C(M 7 ) 
n+1 

M| n + 1, n + 1 2 2\/A „ xa , , „ x9/3 , 2-\/A 



< 



< 



24 _( -~2v A ' 24 v A 
1 

~48 
1 

~96 



(n +!)<«( 



2 1 / 3 (n+l) 2 / 3 (^-l) 
n + 1 

2 1 / 3 (n + l) 2 / 3 (^-l) 
n + 1 



n + 1 

-i 3 



l) 3 + C # (n + l) 2 / 3 ( 



n + 1 



1) + C(M 7 ) 



+ C # (n + l) 2 / 3 ( 



2VA 
n + 1 



1) + C(M 7 ) 



1 3 



■ C(M 7 ). 



The first line follows from the inequality —7 + log 7 + 1 < 



(7-1) 2 



for 1 < 7 < 2 



( note from ( 6.42 ) that C# is independent of Mr). In the second line, we use 
the monotonicity of x(^^- — l) 2 and of (2^/\ — x)~ x in the region n + 1 < x < 
2VA — ^r^-(n + 1) 1 / 3 . In the succeeding lines, we have used the changes of variables 
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in 



x = 2^/Xy, 1 — y = z and 2^/Xz 3 / 2 — s. For the last line, note that 2 1 / 3 (n 
!) 2/3 (^T - l ) ^ M 7> and we require 



(7.8) 



M 7 > V96C*. 



Remark : In estimating the sum in the second line of (7.7) by an integral, the 
monotonicity of the integrand plays a crucial role : we cannot, for example, use an 
estimate of the form (7.5), as the derivative is not sufficiently small. 

For (2), we have < 1 + 2 1 /Hk+i) 2 ' i ' Calculations similar to the previous 
cases (i) , (ii) and (iii) show that 

C(M 7 ) 



/oo 
-2M 7 



(y)dy\ 



+ Ce- Ml < C(M 7 ) 



)l/3 



7 ■ 



The result follows by splitting the sum Ylm ^ ne following regions 2\/A — 
+ < k + 1 < 2VA+^(n+l) 1 /3, 2 ^A+ %fr(n + l) 1 / 3 < fc+K 3^A 
and 3\/A < k + 1 : we leave the detail to the reader. 
Therefore, for 1 1 Mr 



2 i/3 (n +i)2/3 < frf < 1 + ^6, we have 



(7.9) l og 0„(A)<-^ 



2 i/3 (n + 1) 2/3 ( M_i; 



C(M 7 ) 



For case (v), we use the estimation of [ Jol| given in the Lemma 7T (v) below. 

Summarizing, we have 

Lemma 7.1. Let < 65, Sq, 8 7 < 1 and M5, Mq, M 7 > be fixed number. Suppose 
that Sq, M5, Mq and M 7 satisfy conditions (a),(b),(c) and (d) given at the beginning 
of this Section, respectively. Set 

(7.10) < = 2 1 /3 (?1 + i)V3 (l-M.) so that M. = i * 
v/ v / n + 1 2 x /3( n + 1)2/3 

FFe have the following estimates for the large n behavior of <j) n {X) : 

|log</>„(A)| < Cexp(-cn), 
c which n 

M 5 

2 - n+1 - x 2!/3(„ + i)2/3 7 



1. J/0 < < 1 _ 5 

j — n+1 — ° 7 



for some constants C, c which may depend on 5$ 

2 . If 1 < 2j£ < 1 Ms_ 



|log</»„(A)| <Cexp(-ci 3 / 2 ), 

for constants C, c independent of AJ 5 . 
3- If 1 - tkt^ ^273 < < 1 + 1/S/ M6 , 2 /3 , so ffeaf -M e <t< M 6 , there is a 

J 2 1 /3(„ + i)2/3 — „+i — 2 1 /3(n+l) 2/3 ' °' 

constant C(M§) which depends on Mq, and a constant C which is independent 
of Mq, such that 

f2i<,ll(y)dy\ <^^ + Ce^<'\ 
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4. J/1 



Mr 



< 



^<1 + ^ 



2 i/3( n+ i)273 - n+ 



log0„(A) < _^ + C 7(M 7 ), 
96 



for a constant C(Mj). 

+1 



5. ffo§ i/i + ( y 7 < jbg, 



^n(A) < Ce-° A < Ce _cn . 

The really new resul ts in this Lemma are (iii) and (iv). Ind eed, (i) and (ii) can 
also be obtained from (LIS ), and as indicated, (v) is given in [[jolf. 



8. De-Poissonization Lemmas 

In this section, we present two Lemmas which show that 4> n (N) is a good ap- 
proximation Of Qrt.N — fN,n/Nl. 

We need a Lemma show ing the monotonicity of q n ,N in N . The statement and 
proof can be found in [ Jol | . 

Lemma 8.1. For all n,N > 1, 

q n ,N+i < ?n,iV- 

Using this monotonicity result, the following Tauberian-like "de-Poissonization" 
Lemma can be proved. This is a modification of Lemma 2.5 in [Jol] and the proof 
is the same. 

Lemma 8.2. Let m > be a fixed real number. Set — N + (2\/m + 1 + 

l)y/NlogN and = N - (2-y/m + 1 + l)y/N \ogN. Then there are constants 
C = C(m) and Nq — AT (m) such that 

M^) ~ < 1n,N < M^N n) ) + 

for N > N , < n < N. 

The reader will observe that the above Lemma is actually enough for all of our 
future calculations. Nevertheless, for convenience and the purpose of illustration, 
we use the following Lemma for the convergence of moments. 

Lemma 8.3. There exists C > such that 

q n ,N < C(t> n {N - VN), 1 - q n , N < £7(1 - f> n {N + y/N)) 
for all sufficient large N , < n < N . 
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Proof. Note that q n ,N > 0. Using Lemma B.l and Stirling's formula for sufficiently 
large N, we have from fll.llD , 

4> n (N-VN)= 22 7777T1 —Qn,N> 



N'=0 

N 



{N>) 



e -(N-^N) (N _^ )N > 

^ 2^ {m\ qn > N ' 

N>>N--/N 

- qn > N 2^ nvo! 

n'>n-Vn 

>C?n,iV 2^ (jV / ) ^ + l/2 e -7V' =Cqn,N }^ R ' 

N'>N-^Tn N'>N-VIf 

where /(x) = — (N — y/~N) + x\og(N — y/~N) + x — [x + |) logx. One can easily 
check that /(x) is a decreasing function for x > (N — V~N). Thus 

for sufficiently large N, < n < N. 

For the second inequality, note that g„jv < 1 by definition. Again, using 



Lemma 8.1 and Stirling's formula for sufficiently large N, 



1 - 0„(7V + VN) = £ J (1 " ?n,ivO 

>C(1 £ e^') 

where = — + %/iV) + xlog(iV + \/JV) +x— (x+ \) logx. One can check that 
for N < x < N + v/JV, ,g"(x) < so that min^x) = mm(g(N), g(N + VW)). If 
N is sufficiently large, min(g(N), g(N + a/F)) = g(N + VN) = -|log(JV+ 
Therefore 

1 — (f) n (N + V~N) > C(l — g n ,iv )v / iVe 9(A ' ) > C(l - 
for sufficiently large N, < n < N. □ 



9. Proofs of main Theorems 

In this Section, we prove the main Theorems. 



Proof of Theorem 1.1 Assume for defmiteness that i < 0. For t > 0, the 

flf.r 



calculation is similar. From the definition of q n ^ = %T > 



= Prob & J// <t) = q [2VW+tN1/%N . 



N 

Set 



n 



[2VN + tN 1/6 } 
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As t is fixed, observe that < n < N, as N — > oo. Using Lemma 8.2 with any fixed 
value of m > 0, we have 



c 

TV™ 



Set 



o, /,,(™) 

t JV = 2 1 /3 (n + l)2/3 (l _fV^ ) , 

n + 1 



(cf. the definition of f in Q7.1C| ).) Then, for all large N, 

It < tw < — i, and lim tjv = f. 

Let Mq > 2 \t\ be any suffici entl y large, fixed number satisfying condition (c) in 
Lemma [7.l[ Using Lemma 7T (iii), we have, for some constant C(M@) which 
depends on M§, and a constant C which is independent of Mq, 

Mv ( N n) ) = ex p(/°° *™i&(v)dv) (1 + °M 6 (^) + 0(e" W a )) 
Taking — * oo, and then taking Mg — > oo, we obtain, 

J^^MmaT 5 ) = exp (/ 2im i"2(y) d y 

For (finfajlp ), we obtain the same limit by a similar calculation, 



A' 



lim ^(Vjv^) 



exp 



2im{lj 2 {y)dy 



Thus, recalling -p 2i (mf 11 ) 22(2) = w 2 (a;) in (2.18), integration by parts yields 



Um Prob { l JL <t 



N- 



AT1/6 



exp 



2im?li(y)dy ) = F(t) 



□ 



Proof of Theorem 1.2 Integrating by parts, 



Ejv(xw)=/ t m dF N (t) = - mt m - 1 F N {t)dt+ mt™- 1 ^- F N (t))dt 



where Fjv(t) = Prob ( 1n n ^ < t J as in Theorem [□] From Theorem 



LI 



have pointwise convergence of Fjv(t) to F(t). We need uniform control of Fjv for 
large AT. Let M > be a sufficiently large, fixed number and < 5 < \ be a fixed, 
sufficiently small number. 

Set n = [2-//V + tN 1 / 6 }. First consider the case when t < -M. If t < -2N 1 / 3 , 
then Fffjt) = Prob(l N < 2^/N + tN 1 ^) < Prob(l N < 0) = 0. For -2AT 1 / 3 < t < 
—M, (O) and Lemma E| yield 



(9.2) F N (t) = q rhN < C<f> n (N - VN) 
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If -27V 1 / 3 < t < -2(S7V 1 / 3 , when N is sufficiently large, 

2y^V~V¥ > 2VN(1 ~ -^)^ > 2VN(1 - J) j 
n + 1 ~ 2\/]V + tiVVe + 1 ~ 2(1 - S) VN + 1 ~ 2 ' 
Thus, using Lemma 7.1 (v), for large N, 

(9.3) 4> n {N - VN) < Ce- cN < Ce ct \ 

If -26N 1 / 3 < t < -M, 

= 2ViV 



< 



<l + 26 



n+1 2VN-2SVN 
and, using the monotonicity of (2vA — x)/x 1 ^ 3 as a function of x < 2\/\, 

il/3 {2y/N-y/N- (n + 1) 



2 i/s 



> 2 1 / 3 



(n + l)V3 

2 v / iV(l - 4??) 1/2 - (2-//V - MiV 1 / 6 + r 



{2y/N - MN 1 ^ + 1)V3 
as iV -> 00. Thus, we have for -28 N 1 ^ < t < —A/, 



> 



M 



2!/3( n + l) 2 / 3 



n+1 



Therefore, from Lemma 7.1 (iv), provided M- satisfies condition (d) and 28 satisfies 
condition (a), 



2y/N-VN _ ^3 



(9.4) ^(7V-V7V)<CMexpf--i(2 1 / 3 (n+l) 2 / 3 ( 

\ 96 rt + 1 

On the other hand, using the monotonicity of (2y/\ — x)/x 1 ^ 3 as a function of 
x < 2a/A, 

1/3 / 2ViV-VlV-(n + l) \ 

V (n+1) 1 / 3 y 

> 2l/3 f 2Vg(l - ^) 1/2 - (2V^V + ^ 6 + 1) ^ > 



(2VF + tJW« + 1)1/3 



for all -28N 1 / 3 < t < —M, as N -> 00. Therefore (Q gives us 
(9.5) " ~ " - ' 



0„(iV-v / iV)<C(Af)exp(— i 3 ) 



for -2<5iV 1 / 3 < t < -M. 

Inserting the above estimates ( fO| ) and (3.5) into ( |9.2[ ), we obtain 

(9.6) F N {t) < C(M)e c{i for -2iV 1/3 < t < -A/, 

and as Fjy(t) = for £ < — 27V 1 / 3 , it follows by the dominated convergence theorem 
that 



(9.7) 



lim 

N— >oo 



mt m - 1 F N (t)dt = / mt m - l F(t)dt 
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Now consider the case when t > M. If t > iV 5 / 6 - 27V 1 / 3 , then 1 - F N (t) = 
l-Prob(l N < 2VN + tN 1 / 6 ) < l-Prob(l N > N) = 0. For M <t< iV 5 / 6 -2iV 1 / 3 , 
again, ( |9.l[ ) and Lemma |S.3| yield 

(9.8) 1 - F N (t) = 1 - q n , N < C(l - <j> n {N + VN)). 

If 2SN 1 / 3 < t < TV 5 / 6 - 27V 1 / 3 , when N is sufficiently large, 



2VN + VN 2VN(1 



n + 1 



< 



= < 1- ^« 

2VN + 25V N 2 



Thus, using Lemma 7.1 (i) 



(9.9) 1 - cf> n (N + VN) < Ce- cn < Ce~ c ^ < Ce- ctV \ 

If M <t< 25A 1 / 3 , similar calculations to the case -25A 1 / 3 < t < -M yield 



1 r 2Vn + Vn 

- < 1 - 26 < — — — < 1 - - 

2 n + 1 2 1 / 3 (n + l) 2/3 



Therefore, from Lemma 7.1 (ii), as N — ► oo, 



(9.10) 1 - UN ~VN)< expf-c^n + l)2/ 3 ( l - 2y/ ^f N )f^ , 
provided \M satisfies condition (d). However, as in the case -257V 1 / 3 < t < —M, 



we have 



1/a f(n+l)-2y/N + VN 



2 1 / 



> 2 1 / 3 



(2VN + tN 1 / 6 )-2VN(l + ^ w ) 1/2 \ t 



> -. 

2 



(n + 1) 1 / 3 )'-" V (2^N + tN 1 / 6 ) 1 / 3 

for all M <t< 25N 1 / 3 , as N -> 00. Therefore fl9.1Cp gives us 

(9.11) 1 -MN-VN) < Cexp(-ct 3/2 ) 

for M < t < 257V 1 / 3 . 

Inserting the above estimates ( |9.9[ ) and ( p.ll| ) into (9.8), we obtain for M <t < 
TV 5 / 6 - 2JV 1 / 3 

(9.12) l-Fjv(t) < Ce~ c * 3/5 

as TV -> 00. Once again, as 1 - F N (t) = for i > iV 5 / 6 - 27V 1 / 3 , it follows by the 
dominated convergence theorem that 

/>oo />oo 

(9.13) lim / mt m - 1 {l- F N {t))dt= / mi" 1 " 1 ^ - F{t))dt. 

W^oo J J 

□ 



Appendix A. 

As advertised in the Introduction, in this Appendix we give a new derivation of 
the formula 



(A.l) 



E 

N=0 



\ N F N (n) 



det(tZj_fc)o<j,fc< n _i, 
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where dj = L v exp(2-\/Acos6' — ijff)d6, and F^(n) is the distribution func- 

tion for the length, £n(tt), of the longest increasing subsequence in the random 
permutation tt from S^- We set F o (0) = 1. 

Let [i = (/ii, /i2j • ■ ■ , Mr, 0, 0, . . .), fjti > (12 > . . . , be a partition of N, i.e. [ij-, 
1 < j < r, are positive integers and N = /ii + • • • + /i r ; we write /x h N. With /x 
we can associate a Young diagram, also denoted by /x, in the standard way, see for 
example [Sa] . In the Young diagram there are \ij boxes in the j : th row. If we 
insert the numbers 1, . . . , N in the boxes in such a way that the numbers in every 
row and column are increasing we get a (standard) Young tableau t ; t has shape 
/x, s(t) = /i. Let r(ix) denote the number of rows in /x. 

Schensted, [Sc], has constructed a certain bijection, the Schensted correspon- 
dence, between the permutation group Sn and pairs of Young tableaux (t, t') with 
the same shape s(t) = s(t') = [i, where /x h N. This correspondence has the prop- 
erty that if Sn ^ — * (t, t'), fi = s(t), then £/v(7r) equals the length, /ii, of the first 
row in /x, and the length, £^(71"), of the longest decreasing subsequence in ir equals 
r(/i), the number of rows in /x. For details see [Sa]. 

If we put the uniform probability distribution on Sn then clearly the random 
variables In and £' N have the same distribution (just reverse the permutation). Let 
/(/x) denote the number of Young tableaux with shape xx. Then, by the Schensted 
correspondence, 



(A.2) 

If we set hj 
and Young, 

(A.3) 



see for example [Si] 
/ij-i - Hj + 1 > 1. 

(A.4) 



F N (n) = ij 



E 

fj.\-N 
r(/j.)<r, 



m 2 



-j, r = r(/i), we have the following formula, due to Frobenius 



/(/x)=iv! n 

l<i<j<r 

. Note that N = Y%=i Mj : 
Combining the formulas 



«)il^. 



Fjv(n) = Y! 



,-!-./ r(r-l)/2and/i i _i 
) and (Y3) we get 



1 



ipwrip^ 

(*) j=i v J ; 



where the (*) means that we sum over all different integers hi > 1 such that 
Y^hj = N+r(r—l)/2, and A(/i) = Yl i< j(hj — h i ) is the Vandermonde determinant. 
That we can remove the ordering of the hj 's in (A.4) follows from symmetry under 



permutation of hi, 
the Poissonization 

(A.5) <f> n (\) = 

where 



, h r . The constraint Y h-j — N + r(r — l)/2 is removed by 



00 n 

X E Jf\ F M = e " A [! + E A-^-^iMA)], 



= 1 



E 



a(/,) 2 n 



We have used the fact that Y^j > 1 + ■ ■ ■ + r = r(r — l)/2 + r and N > r, since 
the /ij 's are different integers. The condition that the hj 's are different can then 



56 



JINHO BAIK, PERCY DEIFT, AND KURT JOHANSSON 



be removed since otherwise A(h) = 0. Observe that H r (X) is a Hankel determinant 
with respect to the discrete measure 

A m 



v{{m}) = 



(to!) 2 



TO G Z+, 



see [Szl], i.e. 



H r {\) = det( V TO J+fe — - w )o<j.fe<,-i. 

m=l v ; 



If 9j 5 i > 0, are any polynomials with degg^ = j and leading coefficient 1, row and 
column operations on the determinant gives 



(A.6) 



H r (X) = det(^ qj(m)q k (m) 



A' 



(to!)' 



■ )o<j. 



k<r-l- 



We now make a particular choice of qj, qj(x) = x(x — 1) . . . (x — (j — 1)), if j > 1 
and qo(x) = 1, so that 



(A.7) 



d 7 ' 

'-^a m = g J (m)a m , to, j > 0. 



The elements in the Hankel determinants can then be written 



(A.. 



Now, 



Yl 1j(m)qk(m) 



m=l 



(to!) 2 



, k dP d k 
a J b k --- —r > 

rlnJ rlhk Z-^ 



da? db k ^ (to!) 2 



a—b— \/A 



E 



a m b r ' 



^ (to!) 2 2tt 



d0 



and hence we can perform the differentiations in (A.S) and get 
5^ % (m)%(m)-^— = \ {]+k)/2 d^ k - 5 j0 5 m 

m=l ^ 

where dj-k = (27r) _1 exp(2-\/A cos # — i(j — k)9)d9. Inserting this identity into 
the formula flA.6| ) yields 

(A.9) H r (X) = \ r{r - 1)/2 {D r - D r _i), r > 1, 

where D r is the Toeplitz determinant det(dj-k)o<j,k<r-i and Do = 1. Hence, using 
the formula (A.J), we get </>„(A) = e~ x D n , which is what we wanted to prove. 

In the remaining part of this appendix we will give a heuristic argument showing 
why we can expect the random variable £n{k) to behave like the largest eigenvalue 
of a random hcrmitian matrix. From our considerations above we sec that 



fi\-N 



By the same computations as above this leads to 

oo 

(A.10) 0„(A) = e- A [l + ]T \- r ( r -^ 2 H r (\; n)}, 



r=l 
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where 



ft.e{l,...,n+r-l}'' J 

Note that H r (X; n) f H r (X) as n — > oo. We can think of 
(A.ll) 

X hj 1 



-!ffr(A) 



2Ei<, logl'ii-'ijT +E,[(log A)/ij+21og(ftj!)] 



as the probability of the configuration ft, G Zj_. This probability has the form of a 
discrete Coulomb gas on Z+ at inverse temperature (3 = 2 confined by an external 
potential. An N x N random hermitian matrix with a probability density of the 
form Z„ l exp(-TrV(M)) has an eigenvalue density 

J_ e -2Y, i<i lo S \x i -x j \- 1 +Y, :i V(x : i) 



with x € K ; Xi, . . . , a; at are the eigenvalues of M. Thus we can think of the hj ;s 
as some kind of "eigenvalues" . 
Let 

P r (X;n) = H r (X;n)/H r (X), 

i.e. P r (X;n) is the probability that the largest "eigenvalue" is < n + r — 1. Then, 
by (|A]9|) and flA~lC| ), 

(A.12) 

OO OO 

^„(A) = e- A [l + ^P r (A;n)(^ r - A-i)] =e- A + ^P r (A;n)( ( /» r (A)~ ( /» r _ 1 (A)). 

r— 1 r— 1 

Now, the essential contribution to the right-hand side of ( |A.12 ) comes from r around 
since otherwise (f> r (X) — r _i(A) is very small. Thus 

4>n{X) « P 2VX (A;n), 

i.e. </>n(A) is like the probability that the largest "eigenvalue" in the discrete 
Coulomb gas (A.ll) is < n + 2\/A. 



References 

[AS] M.Abramowitz and I.A.Stegun, Handbook of Mathematical Functions, Dover Publications, 
New York, (1965). 

[AD] D.Aldous and P.Diaconis, Hammersley's Interacting Particle Process and Longest Increas- 
ing Subseguences, Prob. Th. and Rel. Fields, 103, 199-213, (1995). 

[Bo] A.Borodin, Longest increasing subsequences of random colored permutations, Electron. J. 
Combin., 6 (1), R13, (1999). 

[BB] R.M.Baer and P.Brock, Natural sorting over permutation spaces, Math. Comp., 22, 385- 
410, (1968). 

[BC] R.Beals and R.Coifman, Scattering and inverse scattering for first order systems, Comm. 
Pure Appl. Math., 37, 39-90, (1984). 

[BDJ] J.Baik, P.Deift and K.Johansson, On the distribution of the leng th of the second row of a 
Young diagram under Plancherel measure, preprint, LANL E-print xiath.CO/990111ii . 

[BR] J.Baik and E.Rains, Symmetrized increasing subseguence problems, in preparation. 

[CG] K.Clancey and I.Gohberg, Factorization of Matrix Functions and Singular Integral Oper- 
ators, Birkhauser, (1981). 



58 



JINHO BAIK, PERCY DEIFT, AND KURT JOHANSSON 



[De] P.A.Deift, Integrable Hamiltonian systems. Dynamical systems and probabilistic methods 
in partial differential equations, 103-138, in Lectures in Applied Mathematics, 31, edited 
by P.A.Deift, C.D.Lcvermore and C.E.Wayne, American Mathematical Society, Providence, 
RI, 1996. 

[DeZel] J.-D.Deuschel and O.Zcitouni, Limiting curves for i.i.d. records, Ann. Probab., 23, 
852-878, (1995). 

[DeZc2] J.-D.Dcuschel and O.Zcitouni, On increasing subsequences of i.i.d. samples, preprint, 
(1997). 

[DIZ] P.A.Deift, A.R.Its and X.Zhou, Long-time Asymptotics for Integrable Nonlinear Wave 
Equations, in Important Development in Soliton Theory, 2nd Edition, edited by A.S.Fokas 
and V.E.Zakharov, Springer- Verlag, to be published. 

[DKM] P.A.Deift, T.Kriccherbauer and K.T-R McLaughlin, New Results on the Equilibrium 
Measure for Logarithmic potentials in the Presence of an External Field, J. Approx. Theory, 
95, no.3, 388-475, (1998). 

[DKMVZ1] P.A.Deift, T.Kriecherbauer, K.T-R McLaughlin, S.Vcnakidcs and X.Zhou, Asymp- 
totics for Polynomials Orthogonal with respect to Varying Exponential Weights, Internat. 
Math. Res. Notices, no. 16, 759-782, (1997). 

[DKMVZ2] P.A.Deift, T.Kriecherbauer, K.T-R McLaughlin, S.Venakides and X.Zhou, Strong 
Asymptotics for Orthogonal Polynomials with respect to Varying Exponential Weights via 
Riemann- Hilbert Techniques, To appear in Comm. Pure. Appl. Math. 

[DKMVZ3] P.A.Deift, T.Kriecherbauer, K.T-R McLaughlin, S.Venakides and X.Zhou, Uniform 
Asymptotics for Polynomials Orthogonal with respect to Varying Exponential Weights and 
Applications to Universality Questions in Random Matrix Theory, To appear in Comm. Pure. 
Appl. Math. 

[DS] P.Diaconis and M.Shahshahani, On the Eigenvalues of Random matrices, J. Appl. Prob. 
31, 49-61, (1994). 

[DVZ1] P.A.Deift, S.Venakides and X.Zhou, The collisionless shock region for the long-time 
behavior of solutions of the KdV equation, Comm. Pure Appl. Math., 47 no. 2, 199-206, 
(1994). 

[DVZ2] P.A.Deift, S.Venakides and X.Zhou, New Results in Small Dispersion KdV by an Ex- 
tension of the Steepest Descent Method for Riemann- Hilbert Problems, Internat. Math. Res. 
Notices, no. 6, 285-299, (1997). 

[DZ1] P.A.Deift and X.Zhou, A Steepest Descent Method for Oscillatory Riemman- Hilbert Prob- 
lems; Asymptotics for the MKdV Equation, Ann. Math., 137, 295-368, (1993). 

[DZ2] P.A.Deift and X.Zhou, Asymptotics for the Painleve II Equation, Comm. Pure. Appl. 
Math., 48, 277-337, (1995). 

[ES] P.Erdos and G.Szekeres, A combinatorial theorem in geometry, Compositio Math., 2, 463- 
470, (1935). 

[FIK] A.S.Fokas, A.R.Its and V.E.Kitaev, Discrete Painleve equations and their appearance in 

quantum gravity, Comm. Math. Phy, 142, 313-344, (1991). 
[FMZ] A.S.Fokas, U.Mugan and X.Zhou, On the Solvability of Painleve I, III and V, Inverse 

Problems, 8, 757-785, (1992). 
[FN] H.Flaschka and A. Newell, Monodromy and spectrum preserving deformations, I, Comm. 

Math. Phy, 76, no.l, 67-116, (1980). 
[FZ] A.S.Fokas and X.Zhou, On the Solvability of Painleve II and IV, Comm. Math. Phy, 144, 

601-622, (1992). 

[Ge] I.M.Gessel, Symmetric functions and P-recursiveness, J. Combin. Theory. Scr. A, 53, 257 
- 285, (1990). 

[GK] I.Gohberg and N.Krupnik, One- Dimensional Linear Singular Integral Equations vol.1 and 

II, Operator theory, advances and applications ; v. 53-54, Birkhauser Verlag, Basel, 1992 
[GW] D.J. Gross and E.Witten, Possible third-order phase transition in the large N lattice gauge 

theory, Phys. Rcw. D, 21, 446-453, (1980). 
[GWW] I.M.Gessel, J.Weinstein and H.S.Wilf, Lattice walks in 7L d and permutations with no 

long ascending subsequences, Electr. J. Combin., 5(1), (1998). 
[Ha] J.M.Hammersley, A few seedlings of research, Proc. Sixth Berkeley Symp. Math. Statist. 

and Probability, Vol. 1, 345-394, University of California Press, 1972. 
[Hi] M.Hisakado, Unitary matrix models and Painleve III, Modern Phys. Lett. A, 11, no. 38, 

3001-3010, (1996). 



LONGEST INCREASING SUBSEQUENCE 



.1!) 



[HM] S.P.Hastings and J.B.McLeod, A boundary value problem associated with the second 
Painleve transcendent and the Korteweg de Vries equation, Arch. Rational Mech. Anal. 73, 
31-51, (1980). 

[Ka] S.Kamvissis, On the Long Time Behavior of the Doubly Infinite Toda Lattice under Initial 

Data Decaying at Infinity, Comm. Math. Phy., 153, 479-519, (1993). 
[Kn] D.E.Knuth, The art of computer programming, vol. 3 : sorting and searching, 2nd ed., 

Addison Wesley, Reading, Mass., 1973. 
[IN] A.R.Its and V.Yu.Novokshenov, The Isomonodromic Deformation Method in the Theory of 

Painleve Equations, Lecture Notes in Math. 1191, Springer- Verlag, Berlin, Heidelberg, 1986. 
[Jol] K.Johansson, The Longest Increasing Subsequence in a Random Permutation and a Unitary 

Random Matrix Model, Math. Res. Lett., 5, no.1-2, 63-82, (1998). 
[Jo2] K.Johansson, Shape fluctuations and random matrices, LANL E-print math. CO/9903134. 
[Jo3] K.Johansson, Transversal fluctuations for increasing subsequences on the plane, preprint, 

1999. 

[JMU] M.Jimbo, T.Miwa and K.Ueno, Monodromy preserving deformations of linear ordinary 
differential equations with rational coefficients, I. General theory and r-function, Physica D, 
2, no.2, 306-352, (1981). 

[Ki] J.-H. Kim, On the longest increasing subsequence of random permutations - a concentration 
result, J. Comb. Th. A, vol. 76, 148-155, (1996). 

[LS] B.F.Logan and L.A.Shepp, A variational problem for random Young tableaux, Advances in 
Math., 26, 206-222, (1977). 

[Me] M.L.Mehta, Random Matrices, Second Edition, Academic Press, San Diago, 1991. 

[Ok] A.Okounkov, Random matrices and random permutations, preprint, 1999. 

[OPWW] A.M.Odlyzko, B.Poonen, H.Widom and H.S.Wilf, On the distribution of longest in- 
creasing subsequences in random permutations, unpublished manuscript. 

[OR] A.M.Odlyzko and E.M.Rains, On longest increasing subsequences in random permutations, 
in preparation. 

[PS] V.Periwal and D.Shevitz, Unitary- Matrix Models as Exactly Solvable String Theories, Phys. 

Rev. Lett., 64, 1326-1329, (1990). 
[Ra] E.M.Rains, Increasing subsequences and the classical groups, Electron. J. of Combinatorics, 

5(1), R12, (1998). 

[Sa] B.Sagan, The Symmetric Group : Representations, Combinatorial Algorithms, and Sym- 
metric Functions, Wadsworth&Books/Cole, Pacific Grove, Calif., 1991. 

[Sc] C.Schcnstcd, Longest increasing and decreasing subsequences, Canad. J. Math., 13, 179 - 
191, (1961). 

[Sel] T.Seppalainen, A microscopic model for the Burgers equation and longest increasing sub- 
sequences, Electron. J. Prob., 1, no. 5, (1996) 

[Se2] T.Seppalainen, Large deviations for increasing sequences on the plane, Probab. Theory 
Related Fields, 112, no.2, 221-244, (1998). 

[Si] B.Simon, Representations of Finite and Compact Groups, Graduate Studies in Mathematics 
vol. 10, American Mathematical Society, 1996. 

[ST] E.B.Saff and V.Totik, Logarithmic Potentials with External Fields, Springer- Verlag, New 
York, 1997. 

[Szl] G.Szego, Orthogonal Polynomials, American Mathematical Society Colloquium Publica- 
tions, Vol 23, 4th Ed, New York, 1975. 

[Sz2] G.Szego, On Certain Hermitian Forms Associated with the Fourier Series of a Positive 
Function, Comm. Seminaire Math de l'Univ. de Lund, tome supplementaire, dedie a Marcel 
Riesz, 228-237, (1952) (or Gabor Szego : Collected Papers - Vol 3 (1945-1972), 270-280, 
Birkhauser, 1982). 

[TW1] C.A.Tracy and H.Widom, Level- Spacing distributions and the Airy kernel, Comm. Math. 
Phys., 159, 151-174, (1994). 

[TW2] C.A.Tracy and H. Widom. Random unitary matrices, permutations and Painleve, 
preprint, LANL E-print tnath.CO/9811154 

[Ul] S.M.Ulam, Monte Carlo calculations in problems of mathematical physics, in Modern Math- 
ematics for the Engineers, E.F.Beckenbach, ed., McGraw-Hill, 261-281, 1961. 

[VK1] A.M.Vershik and S.V.Kerov, Asymptotics of the Plancherel measure of the symmetric 
group and the limiting form of Young tables, Soviet Math. Dokl., 18, 527-531, (1977). 



60 



JINHO BAIK, PERCY DEIFT, AND KURT JOHANSSON 



[VK2] A.M. Vershik and S. V.Kerov, Asymptotic behavior of the maximum and generic dimensions 
of irreducible representations of the symmetric group, Functional Anal. Appl., 19, no.l, 21-31, 
(1985). 

[Wi] H.Widom, personal communication. 

Courant Institute of Mathematical Sciences, New York 
E-mail address: balkScims.nyu.edu 

Courant Institute of Mathematical Sciences, New York 
E-mail address: deiftacims.nyu.edu 

Royal Institute of Technology, Sweden 
E-mail address: kurtjamath.kth.se 



